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Abstract: The problem of improving the rehabihty of perturbative QCD predictions at 
moderate energies is considered. These predictions suffer from substantial renormalization 
scheme dependence, which is illustrated using as an example the QCD effective charge ap- 
pearing in the static interquark potential and the QCD corrections to the Gross-Llewellyn- 
Smith sum rule in the deep inelastic neutrino-nucleon scattering. To reduce this scheme 
dependence, it is proposed to use a modified perturbation expansion, based on a modi- 
fied couplant, which is perturbatively consistent with the conventional running coupling 
parameter, but is free from Landau singularity and has smaller renormalization scheme 
dependence. The modified couplant is obtained by integrating the renormalization group 
equation with an appropriately constructed nonpolynomial generator. The renormalization 
scheme dependence of the perturbative predictions in the modified expansion is discussed in 
detail, including the predictions selected by the Principle of Minimal Sensitivity. It is shown 
that the modified predictions are much more stable with respect to change of the renor- 
malization scheme parameters than the predictions obtained in the conventional approach. 
It is also found that the modified predictions display somewhat weaker energy dependence 
than the predictions obtained with the conventional expansion, which may be interesting 
from the point of view of comparing the low energy and high energy determinations of the 
strong coupling parameter. 
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1. Introduction 



In most phenomenological applications the perturbative QCD predictions are usually eval- 
uated using the renormalization group improved perturbation expansion. This expansion 
works very well at high energies, because of the asymptotic freedom property of the QCD 
running coupling parameter |||, ^. There are however many QCD effects of considerable 
interest, for wich the characteristic energy scale is not very high. Unfortunately, at low and 
moderate energies the predictions obtained from the conventional renormalization group 
improved expansion are not so reliable. This is partially related to the fact that in QCD 
there are many viable definitions of the coupling parameter, corresponding to different 
choices of the renormalization scheme; each choice gives (numerically) a slightly different 
perturbative prediction, as was discussed in ||, |, |, §, 0, % |, |ll|, IH, H, 0, 
In a given order of perturbation expansion this renormalization scheme (RS) dependence 
of perturbative predictions is formally a higher order effect, and its significance diminishes 
at very large energies, where the running coupling parameter becomes small. It turns out, 
however, that in the region of moderate energies — say, of the order of few GeV — the 
renormalization scheme dependence becomes numerically quite large, even with a conserva- 
tive choice of the scheme parameters, as has been pointed out by the author of this article 
in 19, Another problem, which plagues the conventional renormalization group 

improved perturbation expansion at low and moderate energies, is the presence of the so 
called Landau singularity — the running coupling parameter may become infinite at some 
nonzero energy and then below this energy the perturbative predictions simply do not exist. 

Various attempts have been made to improve the reliability of the QCD predictions 
at moderate energies. In the case of physical quantities defined at timelike momenta it 
was observed in [^] that one may reduce the renormalization scheme dependence of the 
predictions by resumming some higher order corrections with the contour integral technique 

|2 



2|,|7|. Other 



in the complex momentum space ||2^, this was verified in detail in 
propositions included a resummation of the series expansion for physical quantities via the 
Pade approximants ||2^, 29, |3^, 31]; modification of the expansion by enforcing certain 

, EoH; and 



analyticity properties in the complex energy plane |32, 33, 0, 35, |3^, 37, 38 



application of more complicated resummation methods |41, 42, 43, 44, 45, 46, 47, 48, 49 
|50| , p2[ . It appears, however, that neither of the proposed methods provides a complete 
and satisfactory solution of the problem. 

The aim of this article is to present an alternative approach. Our starting point is 
the observation that strong RS dependence of perturbative predictions obtained in the 
conventional approach is largely a consequence of the very strong RS dependence of the 
conventional running coupling parameter. We propose therefore to introduce a modified 
running coupling parameter, obtained by integrating the renormalization group equation 
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with a modified, nonpolynomial /^-function. The sequence of modified /3-fiinction is a gen- 
erahzation of the sequence of polynomial /3-function approximants used in the conventional 
expansion. This sequence is constructed in such a way so that the nonpolynomial approxi- 
mants for each order satisfy certain general constraints and ensure the reduced renormaliza- 
tion scheme dependence of the effective coupling parameter. In particular, they are chose 
to ensure the absence of singularity at nonzero positive energy. The modified perturbation 
expansion for physical quantities is then constructed by replacing the conventional cou- 
pling parameter by the modified running coupling parameter. Using as an example the 
QCD effective charge appearing in the static interquark potential and the corrections to 
the Gross-Llewellyn- Smith (GLS) sum rule in deep inelastic neutrino-nucleon scattering, 
we show that perturbative predictions obtained at moderate energies in such a modified 
expansion are much more stable with respect to change of the renormalization scheme. We 
also find that the predictions obtained from the modified expansion have a wealcer depen- 
dence on tlie characteristic energy scale than the predictions obtained form tlie conventional 
expansion. Making some fits to the experimental data for the GLS sum rule we show that 
this effect might have interesting consequences for the QCD phenomenology. 

This article is organized as follows: to prepare the stage for further discussion, we 
briefly summarize in Section 2 tlie main facts about the renormalization scheme dependence 
of the conventional perturbative predictions in QCD. As a concrete example, we discuss 
in detail the RS dependence of the perturbative expression for two quantities, the QCD 
effective charge related to the static interquark potential, and the QCD correction to the 
GLS sum rule, in order to show explicitly that at moderate energies the RS dependence 
is indeed quite substantial. We then briefly discuss previous attempts to improve the 
reliability of perturbative QCD predictions. In Section 3 we describe the construction 
of an improved running coupling parameter that has much weaker RS dependence than 
the conventional coupling parameters. We first give some general constraints that such a 
modified coupling parameter should satisfy in order to give more stable predictions, and then 
we present a concrete model of the modified /3-function. In Section 4 we discuss perturbative 
predictions for physical quantities, obtained with the modified coupling parameter in various 
renormalization schemes. In particular, we consider the defining equations for the so called 
Principle of Minimal Sensitivity (PMS) scheme in the modified expansion, which plays 
important role in our approach. An interesting aspect of our approach is that the sequence 
of the modified /^-functions contains a free parameter, which gives us a natural way of 
utilizing the information of nonperturbative or phenomenological character to improve the 
accuracy of perturbative predictions in low orders of perturbations expansion. We propose 
certain procedure for fixing this parameter, involving phenomenological expression for the 
effective charge related to the interquark potential. In Section 5 we discuss the problem of 
uniqueness of the predictions obtained in the modified expansion: we consider an alternative 
model of the modified running coupling parameter and we evaluate perturbative predictions 
with this coupling parameter. We argue that predictions obtained in the modified expansion 
should be insensitive to the choice of the concrete form of the modified coupling parameter 
even in low orders of the perturbation expansion, provided that the PMS scheme is used. 
In Section 6 we examine, how the use of the modified expansion might affect the QCD 
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phenomenology. We give an argument that for all physical quantities (that belong to the 
general class discussed in this article) the modified predictions obtained in the PMS scheme 
lie below the conventional PMS predictions and evolve less rapidly as a function of the 
characteristic energy scale. Using the QCD correction to the GLS sum rule we illustrate 
that this might be a welcome trend, improving the consistency between low and high energy 
determinations of the strong coupling parameter. In Appendix A we collect some facts 
concerning the renormalization scheme dependence of perturbative QCD predictions; in 
Appendix B we describe the equations for the Principle of Minimal Sensitivity scheme 
in the conventional expansion, using our parameterization of the scheme dependence; in 
Appendix C we give a solution of certain inequality, which is used to select "reasonable" 
scheme parameters in the discussion of the renormalization scheme dependence in Section 
2. 

Basic ideas presented in this article have been previously communicated briefly by the 



author in |53| and 

All the symbolic and numerical calculations reported here have been performed with 
Mathematica. 



2. Renormalization scheme dependence of the perturbative QCD predic- 
tions 

2.1 General form of the perturbative approximants in various renormalization 
schemes 

In order to prepare the stage for further discussion and introduce appropriate notation let us 
briefly recall some basic facts about the renormalization scheme dependence of perturbative 
QCD predictions in finite order of the perturbation expansion. (More details are given in 
the Appendix A.) We shall concentrate on the class of simplest QCD predictions that may 
be expressed in the form of a dimensionless quantity 5 depending on a single variable with 
dimension of (energy)^, which we shall further denote as Q^. We shall also assume that 
the effects of nonzero quark masses are approximated by the step-function Q^-dependence 
of the number n/ of the "active" quark flavors, and we shall restrict our attention to the 
class of mass and gauge parameter independent renormalization schemes. Under these 
assumptions the A^-th order renormalization group improved perturbative expression for 6 
may be written — apart from a multiplicative constant — in the form: 

5W(g2) = a(^)(g2) + ria(^)(Q2) + r2af^^{Q^) + ... + ^N)afN)iQ^)] , (2.1) 

where a(Ar)(Q^) is the A'"-th order effective coupling parameter (the "couplant"), related 
to the gauge coupling parameter g and the strong coupling constant a<j via the relation 
a(iv)(Q^) = 5(7v)(Q^)/4^^ = ai^HQ^)/vr. The couplant 

«(Af)(Q^) satisfies the A^-th order 

renormalization group (RG) equation: 
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Under the assumptions listed above the coefficients b and in the /3-function are ordinary 
numbers, and the coefficients are independent of Q^, so the whole Q^-dependence of 5 
comes then from the Q^-dependence of the couplant a[N){Q'^)- 

In Equation (|2.lD we assumed that the expansion for 6 starts with the couplant a in the 
first power, which is the most common case. The expansion for S may of course begin with 
aP, where p 7^ 1; it is straightforward to generalize our discussion to include this case, but 
we may also note that our considerations could be directly applied to In the following 
we shall usually omit the index in a(^), assuming that 6^^^ is always evaluated with the 
couplant satisfying the A^-th order RG equation. 

The results of perturbative QCD calculations are usually expressed in the modified min- 
imal subtraction (MS) renormalization scheme jH, but there exists a continuum of other 
renormalization schemes, which correspond to different choices of the finite parts of the 
renormalization constants. If we change the renormalization scheme, then in general the 
expansion coefficients rj would change, too. In the new scheme also the coefficients in the 
/3-function may be different. Under our assumptions the coefficients b and ci are indepen- 
dent of the renormalization scheme, but the coefficients Cj for i > 2 are in general scheme 
dependent. The change in the coefficients compensates for the finite renormalization of 
the coupling parameter, but in finite order of the perturbation expansion such a compen- 
sation may be of course only approximate, so the actual numerical value obtained for J^^-* 
in any given order does depend on the choice of the renormalization scheme. The differ- 
ences between values of 5^^^ calculated in various renormalization schemes are formally of 
the order a^+^((5^), but numerically for of the order of few GeV^ they may become 
quite significant, which creates a practical problem when we want to confront theoretical 
predictions with the experimental data. This is the problem that we want to address in 
this article. 

In order to study the renormalization scheme dependence of the perturbative predic- 
tions we need a convenient parameterization of the available degrees of freedom in choosing 
the perturbative approximants. The freedom of choice of the RS in the A^-th order of the 
expansion may be parameterized (in the class of mass and gauge independent schemes) 
by A'^ real parameters. As was discussed for example in ||l9|] (see also Appendix A), it is 
convenient to choose as one of these free parameters the coefficient ri; it is the only free 
parameter in the NL order. For N > 2, following |12|, we shall use as the additional free 



parameters the coefficients C2,---cn in the /3-function. The coefficients r2,---rN are then 
determined from the RS-invariant combinations P2,---PN of the expansion coefficients 
and Ci |11, m, 14, H], beginning with 



P2 = C2 + r2- ciri - rj. (2.3) 

The couplant a(7v)(Q^) in each scheme is determined by the implicit equation of the 
form: 

^1"^^^ = rf^^^^ -n -hciln^ H — + cilnarN) + F'^'^\arN), C2, ■-, cn), (2.4) 



MS 
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where 



C2, 



-Cat) 



da 



2/3 W(c 



a'2 a' 



(2.5) 



This equation is obtained by integrating the renormahzation group equation ( |2.2D with an 
appropriate boundary condition. (See Appendix A for more details.) The role of subtrac- 
tions introduced in the definition of the function F^^^ is to make the integrand finite in 
the limit a — > 0, which simplifies both analytic and numerical evaluation of the integral in 
various cases. The presence of Aj^ and in this general expression may at first seem 
surprising, but it results from the fact that we have explicitly taken into account the exact 
relation (|A.10 ) between A parameters in various schemes, thus choosing Aj^ as a reference 
phenomenological parameter. 

2.2 Propositions for the optimized choice of the renormalization scheme 

In order to resolve the RS ambiguity in finite order perturbative predictions, several methods 
were proposed for making an "optimal" choice of the RS, which we shall now briefly review. 

Conceptually the simplest method was to define the renormalized coupling constant in 
such a way, so as to absorb the radiative corrections to some vertex functions at certain 
configurations of the momenta; in this way one obtains a class of momentum subtraction 
schemes |§, |5|, |6|, 0]. Unfortunately, the vertex functions are in general gauge dependent, so 
instead of the renormalization scheme ambiguity one obtains a gauge parameter ambiguity. 

Another proposition directly referring to diagrammatic calculations and motivated by 



analogy with QED was formulated in |16|: it was proposed to choose the renormalization 
scale in the NL order expression in such a way that the contribution to the physical quantity 
from the vacuum polarization effects due to quarks (represented in the coefficient rj''^'^^ in 
the formula ( |A.5| ) by the term depending on the number n/ of "active" quarks) are absorbed 
into the definition of the renormalized coupling constant (via the n/-dependent term in the 
coefficient h). Extensions of this so called BLM procedure to the case of scale fixing in 
the presence of some higher order corrections were discussed in 



56, 57 



5G 



Unfortunately, it proved difficult to extend the BLM method to higher orders in such a way 
that it would lead to a unique choice of all the renormalization scheme parameters relevant 



in the given order of the perturbation expansion, as was discussed in ||6l|, |6^, |65|| . 

A more radical proposition, originating from the work of |^, was to choose scheme 
(by adjusting the constants in the Equation ( [A.7|) ) in such a way that all the expansion 
coefficients in the A'^-th order approximant 5^^^ for the physical quantity 5 are identically 
zero, r-i = 0, without any reference to any explicit subtraction procedure. In other words, 
in this scheme the renormalized couplant coincides with the physical quantity 8. In our 
parameterization this corresponds to the choice of the scheme parameters r\ = 0, Ck = Pk- 
This choice was dubbed in [11| the Fastest Apparent Convergence (FAC) scheme. 

Another interesting proposition was to choose the scheme parameters according to the 
so called Principle of Minimal Sensitivity fl^, 11, |l^ (for a more recent discussion see ||66| ): 



since the physical predictions of the theory should in principle be independent of the choice 
of the renormalization scheme, we should give preference to those values of the scheme 
parameters, for which the finite order perturbative predictions are least sensitive to the 
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local changes in these parameters. In the NL order the scheme parameter ri corresponding 
to the PMS scheme is a solution of the equation: 

^(5«(a(Q2,ri),ri) \r^=f,= 0, (2.6) 

where we emphasized the fact that 5^^^ depends on ri both explicitly and implicitly via the 
ri-dependence of a(Q^). In the NNL order the parameters fi and C2 corresponding to the 
PMS scheme are solutions of the system of two equations: 

d 

— 5(2)(a((3^ri,C2),ri,C2) |ri=fi,c2=c2= 0, (2.7) 



5(^)(a((5^,ri,C2),ri,C2) |ri=fi,c2=S2= 0, (2.8) 



and 

_d_ 

dc2 

where it is understood that in performing partial differentiation of 5^"^^ over ri and C2 also 
of the implicit dependence of a((5^,n,C2) on these parameters is taken into account. A 
more explicit form of the NL and NNL equations for the PMS parameters is given in the 
Appendix B (it should be noted, that the parameterization of the RS dependence and the 
form of the NNL order scheme invariatn p2 adopted in this article are different from those 



assumed in the original paper |11|). 

The choice of the PMS scheme has nice conceptual motivation, but it may be also 
advantageous from the point of view of resummation of the perturbation series. As is well 
known, the perturbation series in QCD is only asymptotic (i.e. its radius of convergence 
is equal to zero) and a straightforward summation of successive corrections in any fixed 
renormalization scheme must inevitably give infinite result. One may nevertheless give 
meaning to the sum of such divergent series by applying various summation methods, such 
as the method based on the Borel transform. (The problem of large order behaviour in 
various field theory models and appropriate summation methods has been reviewed for 
example in |^ |6§].) In it was conjectured that the PMS method — which involves a 
"floating" choice of the scheme parameters, depending on the energy and the order of the 
expansion — could provide an "automatic" resummation method, which would convert a 
divergent series expansion into a convergent sequence of nonpolynomial approximants. This 
idea was further discussed in ||7^, 73, 74]. Although in the case of realistic QCD 

series this conjecture is far from being proved or disproved, it has been known for a long 
time that re-expansion procedures based on the introduction of some auxiliary parameters, 
which are "floating", i.e. they are fixed in an order dependent way, could indeed provide an 
effective resummation method for divergent series. For example in |75] a rigorous proof has 



been given that a method based on an order dependent mapping gives convergent results in 
the case of the factorially divergent expansion of a one dimensional non-Gaussian integral, 
and a compelling numerical evidence has been obtained that the same is true in the case 
of the divergent perturbation expansion for the ground state of the anharmonic oscillator; 
in a numerical study of the anharmonic oscillator presented in [[T^ ] it has been shown, that 
a sequence of approximants constructed according to the minimal sensitivity criteria seems 
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to be convergent; finally, in |77, |78| a method for improving the perturbation series for the 
anharmonic oscillator has been proposed, which involves order dependent choice of certain 
parameters, and it has been proved, that it results in a convergent sequence. For all these 
model divergent series rigorous proofs have been obtained, that the resummation methods 



based on the PMS criteria do indeed give convergent results 81, |8^. (A more 

general approach to the problem of resumming the divergent series with a nonpolynomial 
sequence of approximants with auxiliary parameters has been presented in |^, |8^, |5|, |86| . 
The author is grateful to Prof. Yukalov for bringing these references to his attention.) 

Regardless of the ingenuity of various propositions for the optimal choice of the renor- 
malization scheme one cannot escape the fact that they are (at least so far) of heuristic 
character. In the scheme parameter space in the vicinity of the "optimal" schemes we have 
many other schemes, which a priori look equally reasonable; predictions in such schemes 
also should be somehow taken into account. Unfortunately, even if we restrict ourselves to 
"reasonably" looking schemes, the variation of the predictions obtained from the conven- 
tional perturbation expansion at moderate energies appears to be quite large. We illustrate 
this problem in the following subsections. 

Before we move on, let us comment on a completely different strategy of evaluating 
perturbative QCD predictions, i.e. the so called method of effective charges, originating 
from the works of [^, ^ U, [l^ and recently further developed in |87, In its simplest 
form this approach is based on the manifestly RS independent evolution equation for (^(Q^), 



^d5m _ _K2 



k=2 



(2.9) 



which is obtained by taking the limit — > 00 in the equation (|A.12|) . The great advantage 
of this approach is that the generator of the evolution for d{Q^) is a scheme independent 
object, so all the scheme dependence of perturbative predictions is apparently absent in 
this formulation from the very beginning. Let us note, however, that if we use as an 
expression for the generator a simple truncated series in the given order, we end up with 
the expression coinciding with the formula obtained in the FAC scheme; this expression 
cannot be considered completely satisfactory, for example because it suffers from the Landau 
singularity problem (at least for physical quantities with positive pk). On the other hand, if 
we try to improve the series expansion for the generator via some sequence of nonpolynomial 
approximants, then in the low orders we encounter the problem of arbitrariness in the 
choice of these approximants, which has a similar effect on the predictions as the effect of 
arbitrariness in the choice of the renormalization scheme in the usual expansion in terms 
of the effective couplant. Another drawback of this approach is that there are physical 
quantities, which cannot be easily expressed in terms of simple effective charges. It seems, 
therefore, that this approach cannot compete at present with the commonly used expansion 
in terms of an effective couplant. 

Another interesting manifestly scheme independent approach has been developed in 



| 89| , |9Q| , |91[| , where authors propose to make a direct comparison of physical quantities. The 



attractive feature of this approach is that by expanding one physical quantity in terms of 
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another physical quantity one may avoid to some extent the compHcations arising from the 
scheme dependence. It was also observed that relating physical quantities at appropriately 
chosen energy scales — called by authors "commensurate scales" — one may eliminate large 
contributions to the expansion coefficients that arise from the terms that explicitly depend 
on the /3-function coefficients. Unfortunately, results obtained with commensurate scale 
relations cannot be easily combined with the results obtained in the more conventional 
approaches to the QCD phenomenology. In particular, in order to achieve their goals. 



authors of ||89| adopt a somewhat unusual multi-scale method, evaluating each term of the 



series expansion at a different energy scale. 



2.3 Renormalization scheme dependence of 5y 



The fact that finite order perturbative predictions obtained with the conventional pertur- 
bation expansion exhibit at moderate energies a strong dependence on the choice of the 
renormalization scheme, despite a conservative choice of the scheme parameters, has been 



demonstrated in several articles ||17|, |18|, [19|, ^ |29||. As a preparation for further 



discussion we illustrate this effect once again, using as an example the perturbative QCD 
expression for the effective charge 5y, appearing in the static interquark interaction poten- 



tial [92|. To our knowledge, a complete analysis of the RS dependence of 5\ has not been 
performed so far. This effective charge is very important for the study of heavy quarko- 
nia, and it would play an important role in our construction of an improved perturbation 
expansion. (The present status of the theory and phenomenology of heavy quarkonia has 
been reviewed in ||93|| .) 

The static interquark potential may be defined in a nonperturbative and gauge invariant 
way via the vacuum expectation of rectangular Wilson loop of size r in spatial dimensions 
and size r along temporal axis |92|: 

V{r) = - lim — In < 0|TrPexp (i5 / dx^^A^T" ) |0 >, (2.10) 

T-*00 IT \ J J 

where V denotes a path ordering prescription. The effective charge dy enters the Fourier 
transform V{Q'^) of V{r) in the following way: 

V{Q') = -A.'c/-^, (2.11) 

where in the SU(N) gauge theory with fermions in the fundamental representation Cp = 
Tf{N'^ — l)/N^ with Tp defined by the normalization of the gauge group generators. 



Tr(r"r'') = TfS"". Perturbative expression for 6y has the form (|2.l|) and is presently 



completely known up to and including the NNL order ||9|, |95|, g^, ^, |99|, |100|, |10]J, and 
some terms in the expansion are known even beyond this order ||102| , tl03| . (Recent theo- 
retical studies of various contributions to Sy have been summarized in ||104l ) For nj = 3, 
which is the most interesting case from the phenomenological point of view, we have in 
the_MS scheme r[°^^^ = 1.75 and r!^^^^ = 16.7998, as well as b = 9/2, ci = 16/9 and 
(MS ^ 4 471^ ^j^ipj^ impHes = 15.0973. 
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In Figure [I] we show <5v as a function of ri , for several values of C2, at = 3 GeV^; the 
NL order prediction is also shown for comparison. (Here and in all other numerical calcu- 

(3) 

lations described in this article we assume — unless stated otherwise — A.)-J- = 0.35 GeV. 

MS 

If we accept the world average 

(M|) = 0.1182 ± 0.0027 obtained in and convert 

this parameter into Al^i., using the procedure described in detail in Section 6, we obtain 

^Ms ~ O-'^^ GeV±Q;Q4.) As we see, the differences between predictions obtained for differ- 
ent values of the scheme parameters ri and C2 are quite substantial. In Figure || we show 
the NNL order predictions for 5\ as a function of Q^, in several renormalization schemes, 
including the MS scheme and the PMS scheme. As we see, for of the order of few GeV^ 
the differences between predictions in various schemes are quite large, although they of 
course rapidly decrease with increasing . 




2.4 The problem of reasonable scheme parameters 

One could of course raise an objection that large differences between predictions obtained in 
different schemes are a result of a wrong choice of the scheme parameters. One could argue 
that if we would choose "reasonable", "natural" scheme parameters, then the differences 
between the schemes would not be very big. In order to verify, whether this indeed might 
be the case we need to give more precise meaning to the intuitive notion of a "reasonable" 
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Figure 2: as a function of (for n/ = 3), as given by the conventional NNL order expansion in 
several renormalization schemes, including the MS scheme (dash-dotted line) and the PMS scheme 
(soHd Hue). 



or "natural" renormalization scheme parameters. Intuitively, a reasonable renormalization 
scheme is a scheme in which the expansion coefficients for the physical quantity b and the 
/3-function are not unnaturally large. It was observed in 19, that one could assess 



the "naturalness" of the renormalization scheme by comparing the expansion coefficients 
in this scheme with their scheme invariant combinations pk, relevant for the considered 
physical quantity in the given order of perturbation theory: a scheme could be considered 
"natural", if the expansion coefficients are such that in the expression for the invariants pk 
we do not have extensive cancellations. The degree of cancellation in pk may be measured 
by introducing a specific function of the scheme parameters, which in the simplest variant 
may be chosen to be the sum of the absolute values of the terms contributing to pk- In the 
NNL order we have: 

o'2{ri,C2) = \c2\ + \r2\ + ci|ri| + rf. (2-12) 

If we choose a scheme which corresponds to the scheme parameters that give rise to extensive 
cancellations between various terms contributing to p2, then the sum of the absolute values 
of the terms contributing to p2 would be much larger than \p2\. The idea is then to estimate 
the magnitude of RS dependence by comparing the predictions evaluated for the schemes 
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that have comparable degree of "naturalness" — i.e. that involve comparable degree of 
cancellation in p2. The problem of selecting natural schemes is in this way reduced to the 
problem of deciding, what degree of cancellation is still acceptable. As was pointed out in 
[l9[| , one may answer this question by referring to the PMS scheme. Using the approximate 
expressions ( [B.14| ) and ( B.15 ) for the parameters in the PMS scheme we find that in the 



weak coupling approximation a2(PMS) « 'AP'A- Therefore, if we accept the PMS scheme 
as a "reasonable" scheme — which seems to be a sensible condition — then we should also 
take into account predictions in the whole set of schemes, for which the scheme parameters 
satisfy the condition 

cj2<2|p2|. (2.13) 

An explicit description of the set of scheme parameters satisfying the condition CJ2 < l\p2\, 
where I is some constant {I > 1) is given in the Appendix C. (Incidentally, this shows, why 
it is important to calculate the NNL order corrections to physical quantities — this is the 
first order, in which we may introduce a constraint on the scheme parameters based on the 
scheme invariant combination!) 

In Figure || we show the contour plot of the NNL order prediction for 5y at = 9 GeV^ 
as a function of the scheme parameters ri and C2. The boundary of the set of the scheme 
parameters satisfying the condition ( ^.13 ) is indicated with a dashed line. One may easily 



verify that scheme parameters for all the curves shown in Figure |^ lie within this region. 
This suggests that strong renormalization scheme dependence of the NNL order perturbative 
predictions for 5\ cannot be simply attributed to the improper, "unnatural" choice of the 
renormalization scheme. 

Let us comment at this point that the frequently used procedure of estimating the 
scheme dependence by varying the renormalization scale parameter /i^ in some range, say 
(5^/4 < /x^ < 4(5^, may give misleading results, as has been pointed out for example in |17|. 
This is explained in more detail in Appendix A. 

2.5 Renormalization scheme dependence of 5gls 

As an another example of a physical quantity of considerable interest at moderate energies 
let us consider the QCD correction to the Gross-Llewellyn- Smith (GLS) sum rule for the 



non-singlet structure function xF'i{x, Q ) in the i/(z^)-nucleon deep inelastic scattering [[1Q6| : 



/' dx \ \f^\x, Q2) + Fl^{x, g2)l =3-3 [<5gls(Q') + ^llsiQ'')] , (2.14) 

where Aq^s is the nonperturbative (higher twist) contribution, and (5gls denotes the per- 
turbative QCD correction. The GLS integral is interesting, because it proved possible 
to measure it directly over large range of Q^, without performing any extrapolation of 
xF3(x,Q^) over . The latest comprehensive analysis of the experimental data for the 
GLS sum rule has been reported in |107|] . The perturbative correction (5gls has the form 
( |A.1|) and is presently known up to the NNL order j^, |108| , |109| , [L10| : for = 3 we have 

in the MS scheme rf^^^ = 43/12 and r\^^^^ = 18.9757, which imphes p^^^ = 4.2361. The 
comparison of the conventional predictions for <5gls in the MS, PMS and EC schemes has 
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Figure 3: Contour plot of the conventional NNL order predictions for Sy at = 9GeV^ (for 
71/ = 3) , as a function of the scheme parameters ri and C2 . The dashed line indicates the boundary 
of the set of parameters satisfying the condition ( 2.13| ) . 



been presented in |111], and the problem of improving the perturbation expansion for this 
quantity has been discussed in | 112 , 113 1. 

We begin our discussion of the RS dependence of (^gls by considering the the contour 
plot of the NNL order prediction for (5gls at = 5 GeV^ as a function of the scheme 
parameters ri and C2, as shown in the Figure ^. The boundary of the set of the scheme 
parameters satisfying the condition ( 2.13 ) is indicated with a dashed line. This set is of 
course smaller than the corresponding set for 6y, because ^2^^^^ is smaller than p^. Let us 
also observe that in the case of (^gls the parameters of the MS scheme fall well outside the 
set singled out by the condition (|2.13| ), since a2{ri^^ , C2^^) /\p2^^\ = 10.07; nevertheless the 
predictions in this scheme do not show any pathological behaviour. 

In Figure ^ we show the dependence of dcLS on the parameter ri at fixed value of 
= 3GeV^, for several values of C2. 

In Figure ^ we show the NNL order predictions for (5gls as a function of (for nj = 3) 
in several renormalization schemes. Since the set of the scheme parameters satisfying the 
condition ( p. 13 ) is in the case of (5gls relatively small, we have included in Figure |6| also the 
predictions for scheme parameters lying slightly outside this set, in order to obtain a better 
picture of the scheme dependence. As we see, also for this physical quantity the differences 
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Figure 4: Contour plot of the conventional NNL order predictions for (5gls at Q2 ^ 5 GeV^ (for 
71/ = 3) , as a function of the scheme parameters ri and C2 . The dashed line indicates the boundary 
of the set of parameters satisfying the condition ( 2.13| ) . 



between the predictions in various schemes satisfying the condition ( ^.13|) become quite 
substantial for of the order of few GeV^. 

It is of some interest to compare the perturbative contribution (5gls at various with 
the estimate of the nonperturbative contribution Aqls(Q^)i which has been discussed in 



several papers [|114| , |115| , |116| . Following |107|] we shall assume: 

HT /^2^ _ (0.05 ±0.05) 

As we see, the differences between predictions in various schemes become large even at 
those values of Q^, for which the nonperturbative contribution is estimated to be small. 



^lUQ') = - GeV^. (2.15) 



2.6 Possible significance of large scheme dependence at moderate 

The discussion of 6\ and 5gls shows that for of the order of few GeV^ the predictions 
obtained with the conventional renormalization group improved perturbative approximants 
are surprisingly sensitive to the choice of the renormalization scheme, even with a conser- 
vative choice of the scheme parameters. The question is now of a proper interpretation of 
this fact. A "pragmatic" point of view — assumed in many phenomenological analyses — 
is to dismiss this fact as some technical "oddity" and continue to use the MS scheme. This 
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Figure 5: Jqls at = 3 GeV^ (for n/ = 3), as a function of ri, for several values of C2, as given 
by the conventional perturbation expansion with = 0.35 GeV. Dashed line indicates the NL 
order prediction. 



is presumably a healthy practical approach, justified to some extent by the experience that 
variation the scheme parameters in the immediate vicinity of the MS parameters does not 
lead to dramatic scheme dependence — at least in those cases, for which the radiative cor- 
rections have been evaluated, because there is no proof that this would be true for arbitrary 
physical quantity. (However, the important question for the applications in phenomenology 
is not only what is the value of the prediction, but also what is the estimated accuracy of 
the prediction — the plots we have shown suggest that perturbative results obtained in the 
MS scheme are less precise than commonly thought.) However, from a theoretical point of 
view such an approach cannot be considered completely satisfactory. 

A different and a rather extreme point of view would be that the observed scheme de- 
pendence is a signal of a genuine breakdown of perturbation expansion, which indicates the 
need for a fully nonperturbative approach already at the moderate energies. Fortunately, 
this seems unlikely to be the case. We see — for example in the case of the GLS sum 
rule — that the energies, at which the differences between predictions in various schemes 
start to become large are higher than the energies, for which — according to the avail- 
able estimates — the essentially nonperturbative corrections may become very important. 
This suggests, that strong renormalization scheme dependence has little to do with the 
unavoidable breakdown of perturbation expansion at very low energies. 



- 15 - 




Figure 6: (5gls as a function of (for n/ = 3), as given by the conventional NNL order expansion 
in several renormalization schemes, including the MS scheme (long-dash-dotted Hue) and the PMS 
scheme (solid Hue) . The estimate of nonperturbative contribution is shown for comparison (short- 
dash-dotted line). 

There is then a third possibihty that the strong renormahzation scheme dependence in- 
dicates the need to abandon the conventional perturbative approach in favor of an improved, 
less scheme dependent formulation, more suitable for moderate energies. One should keep 
in mind that the conventional renormalization group improved expansion is after all just a 
specific resummation of the original, "non-improved" expansion. The substantial RS depen- 
dence of the conventional approximants may simply show that this resummation, despite 
being extremely useful at high energies, is not very well suited for moderate energies. This 
is the point of view that we want to pursue in this article. 

2.7 Previous attempts to modify the renormalization group improved expansion 

The problem of improving the stability of perturbative predictions was considered by many 
authors. It turned out that in the case of physical quantities evaluated at timelike mo- 
menta — such as the QCD corrections to total cross section for the e+e~ annihilation into 
hadrons or the total hadronic decay width of the r lepton — surprisingly good results 
may be obtained by resumming to all orders some of the so called vr^ corrections, arising 



from the analytic continuation from spacelike to timelike region ||2l|, |2^, ^ Such 
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resummation is conveniently performed by expressing these physical quantities as contour 
integrals in the complex-Q^ space, which are then evaluated numerically [|^, 23, 21|. For 
example in the case of the e"'"e~ annihilation into hadrons the ratio R(,+e- of the hadronic 
and leptonic cross sections 



o-tot(,e^e 



hadrons) 



o-tot(e"' 



may be expressed as 



i?g+g_(s) = 127rlmn(s + ie) = — 



da 



D{a) 



c 



a 



(2.16) 



(2.17) 



where n(g^) is the transverse part of the correlator of electromagnetic currents and D{q^) 
is defined as 

D{q^) = -12vr2g2_ln(g2) =^Y.Q)[l + 8u{-q')] • (2.18) 

/ 



{Qf are the charges of the "active" quarks and 5d denotes the perturbative QCD correction, 
which has the form ( [A.l| ).) The contour C runs clockwise from a = s — ie, around a = 
to a = s + ie. 

However, the case of physical quantities evaluated at spacelike momenta — such as the 
QCD corrections to the deep inelastic sum rules — proved to be a more difficult problem, 
home authors |8|, |2|, |0|, 0] advocated the use of Fade approximants to resum the series 
expansion (|2.l|) for the physical quantity 6, observing some reduction of the RS dependence. 
Another group of authors pursued an approach stressing the importance of ensuring the 
correct analyticity properties of perturbative approximants in the complex Q^-plane. The 
basic idea of this so called Analytic Ferturbation Theory (AFT) is to modify the pertur- 
bative predictions by expressing them as dispersion integrals over the timelike momenta, 
with the spectral density chosen in such a way so as to ensure the consistency with the 



weak coupling perturbation expansion |32]. In the leading order this is equivalent to the 
"analytization" of the running coupling constant [|^, 35, |3^, i.e. replacing the lowest order 
expression a(Q^) = 2/(61nQ^/A^) by an "analytic" couplant 



+ 



A' 



lnQ2/A2 A2-Q 



(2.19) 



which is nonsingular for = A^. It was shown that using this approach one obtains 



substantial reduction of the RS dependence in the case of the Bjorken sum rule |37] and the 



GLS sum rule |112]. Within AFT one may also obtain some improvement for quantities 

|. (The AFT approach has been summarized in 



3S 



evaluated at timelike momenta 
p3[| , and it has been extended to a more general type of physical quantities in |117, 118, 119 
|120| ].) The problem of reducing the RS dependence was also discussed within a somewhat 
involved approach formulated in 



111, ^, M- 



Unfortunately, none of these approaches is completely satisfactory. (Besides the efforts 
mentioned above there were also other works, in which authors studied how the appli- 
cation of various resummation methods affects the stability of perturbative predictions 
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1 44, 50, 113, ^ |5^, 45, 47 1, but these attempts concentrated only on the renormaliza- 



tion scale dependence. Unfortunately, such an analysis does not give a proper picture of 
the full renormalization scheme ambiguity of perturbative predictions at higher orders, as 
we tried to explain in the previous subsections.) The Pade approximants, being rational 
functions, tend to develop singularities for some particular values of the scheme parame- 
ters. It also seems that the character of divergence of the QCD perturbation expansion at 
high orders is such that it cannot be resummed by approximants of this type |67, 68]. On 
the other hand, the "analytization" procedure introduces some essentially nonperturbative 
contributions (required to cancel the unphysical singularities of the conventional perturba- 
tive approximants — as may be seen in the formula shown above), which generate l/Q^ 
corrections at moderately large Q^. These corrections are troublesome, because for some 
physical quantities — such as the i?e+e- ratio — they are incompatible with the opera- 
tor product expansion, and for other quantities — like 5gls — such corrections interfere 
with the available estimates of nonperturbative contributions, causing a double counting 
problem. The "analytization" procedure also introduces a very strong modification of the 
low-Q^ behaviour, making all the expressions finite at = — and this modification is 
very "rigid", i.e. it is uniquely defined by the weak coupling properties of the theory, which 
is slightly suspicious. Finally, if we decide to apply the „analytization" procedure directly to 
physical quantities, then we effectively abandon the concept of expressing all the predictions 
of the theory in terms of one universal effective coupling parameter. 

In this note we make an attempt to formulate an alternative method for improving 
perturbative QCD approximants, which gives more reliable results than the conventional 
expansion and yet at the same time does not depart too far from the standard perturbative 
framework. 



3. The modified couplant 

3.1 Strong scheme dependence of the conventional couplant 

A closer look at the conventional RG improved expansion immediately reveals that one of 
the main reasons for the significant RS dependence of finite order predictions is the very 
strong RS dependence of the running coupling parameter itself. An extreme manifestation 
of this scheme dependence is the fact that in a large class of schemes the couplant becomes 
singular at finite nonzero Q^, with the location and the character of the singularity depend- 
ing on the choice of the scheme. In the NL order the running coupling parameter becomes 
singular at = (Qnl)^' where: 



Qnl - ^MS ( ^ 



exp 



ri 



(3.1) 



and for close to (Qnl)^ it behaves as 



a(i)(Q^ 



bci In 



{Q 



NL) 



(3.2) 
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In the NNL order the running coupling parameter has an infrared stable fixed point for 
C2 < 0, while for C2 > it becomes singular at 
the location of singularity is given by: 



(Qnnl) ! fo'^ example, for 4c2 — > 



Qnnl = ^MS 




r\ + Cl In 



+ 



2C2 



2i/c2 4c2 - c- 



: arctan 



Cl 



(3.3) 



and for close to (Qnnl)^ the couplant behaves as 



i{2) 



hc2 In ■ 



(3.4) 



(Qnnl) 

In higher orders we find similar type of behavior, depending on the signs of Cj, i.e. either 
the infrared stable fixed point or a singularity at nonzero (positive) Q^, with the character of 
the singularity determined by the highest order term retained in the /3-function. Although 
the singularity itself occurs in the range of which is normally thought to belong to the 
nonperturbative regime, its strong RS dependence affects also the behaviour of the couplant 
at above the singularity, resulting in a strong RS dependence of the running coupling 
parameter for those values of that lie in the perturbative domain, which cannot be 
properly compensated by the corresponding changes in the expansion coefficients rj. It 
seems unlikely that one would be able to improve stability of the perturbative predictions 
in a significant way without somehow solving this problem. 

Fortunately, within the perturbative approach we have some freedom in defining the 
actual expansion parameter. The idea is then to exploit this freedom and try to construct 
an alternative couplant, which would be much less scheme dependent than the conventional 
couplant, hopefully giving also much less RS dependent predictions for physical quantities. 
From the preceding discussion it is obvious that first of all we would like this modified 
couplant to be free from the Landau singularity. 

3.2 General properties of the modified couplant 

The idea of modifying the effective coupling parameter a{Q'^) in order to remove the Landau 



singularity has of course a long history, dating back to the early days of QCD ||l2l| , |122 



[12^ , m, m, |T2|, |i23, m, 112^, m, im, im, iin, [m, im, m, [m, im, im, mi 
im, m, [m, im, iii^, in^, [ni, m^, 115^, imi, [m, m, im, iin, m, m, m. it tums out, 

however, that various proposed models — although very inspiring — are not directly useful 
from our point of view. Let us formulate more explicitly the conditions that the modified 
couplant should satisfy in order to be suitable for our approach. (The ideas presented here 
have been first formulated briefiy by the present author in p3| .) 

Our first constraints are related to our assumption that we want to stay as close as 
possible to the usual perturbative framework. This means in particular that we want to 
retain the concept of an effective coupling parameter satisfying the renormalization group 
equation, although the conventional, polynomial generator in this equation would 

have to be replaced by an appropriately chosen nonpolynomial function j3'^^\ First of all, 
in order to ensure the perturbative consistency of the modified perturbation expansion with 
the conventional expansion in the N-th order we shall require 
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(I) /3W(a)-/3W(a) = 0(a^+3)_ 

If this condition is satisfied, then the expression obtained by replacing the conventional 
couplant a(Ar)((5^) in the expansion for 5^^^ by the modified couplant a(Ar)((5^) (we may 
temporarily denote such an expression as 6^^^) differs from the original expression only 
by terms that are formally of higher order than the highest order term retained in the 
conventional approximant. 

Secondly, in order to ensure that the modified running coupling parameter remains 
nonsingular for all real positive Q^, we shall require that: 

(Ila) /3(^)(a)isne gative and nonsingular for all real positive a and for a — > +00 it behaves 
like ^(^)(a) ~ -Ca^ where ^ is a positive constant and k < 1, 



or 



(lib) P^^^ (a) has zero for real positive ao and is negative and nonsingular for < a < 



Solving explicitly the RG equation in the case (Ila) with the asymptotic form of the gener- 
ator, we find in the limit ^ for /c < 1: 



A2 

Q- 



1 

i-fe 



a(Q^)~( In -2) , (3.5) 



while for /c = 1 we obtain 



<Q') ~ [q2) ■ (3-6) 

We could have formulated the condition (Ila) in a much stronger way, strictly enforcing 
some type of low-Q^ behavior of the couplant. However, our aim is to improve the QCD 
predictions at moderate energies, where the perturbation expansion is still meaningful and 
gives dominant contribution, although its application may be nontrivial; we do not intend 
to make any claims about the predictions at very low Q^, where noperturbative effects dom- 
inate and a proper treatment requires much more than just an improved running expansion 
parameter. Therefore the exact asymptotic behaviour of the couplant a{Q'^) at very low 
energies is not very important for our considerations; what we basically need is that the 
modified coupling parameter does not become singular at nonzero positive Q^. 

Our next condition is related to the observation that in general it is not difficult to 
obtain a running coupling parameter without the Landau singularity by introducing essen- 
tially nonperturbative (i.e. exponentially small in the limit a ^ 0) terms in the /3-function; 



indeed, several models of this sort have been discussed in the literature |127, |129| , |133| , 34 



35, 36, 151 1 . Unfortunately, such models usually generate corrections at large Q 



2 



which are unwelcome. As we have already mentioned, if such corrections are present, the 
existing estimates of nonperturbative effects based on the operator product expansion be- 
come useless because of the risk of double counting. The l/Q^ correction in the large-Q^ 
expansion is particularly unwelcome, because for many physical quantities there is no room 
for such a term in the operator product expansion. It seems therefore desirable to consider 
only those nonpolynomial functions f3^^^ that satisfy the condition: 
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(III) /3(^)(a) is analytic in some neighborhood of a = 0. 

Analytic /3-functions are preferable also from another point of view: the sequence of per- 
turbative approximants in the modified perturbation expansion with such a generator may 
always be interpreted as a pure resummation of the original expansion, without any terms 
being added by hand. 

Our next constraint is related to the following observation: the finiteness of a{Q'^) for 
all real positive does not automatically guarantee that at moderate (but not very low) 
— corresponding to moderately large a — the couplant would grow less rapidly compared 
to the conventional couplant and that it would be less RS dependent in this range. (It is 
easy to construct explicit counterexamples.) To ensure such a behaviour we shall impose 
an additional constraint on f3^^\a), which may be loosely formulated in the following form: 

(IV) /3(^)(a) - /3(^)(a) = K^^\ci,C2, Civ)a^+^+P + 0{a^+^+P), where p is some pos- 
itive integer and the coefficient K^'^^ is positive for = 1 and for > 2 it is a 
slowly varying function of its parameters, with the property that K^^^ > when the 
coefficients C2,...CAr are large and positive and K^^'^ < for C2,...C7v large negative. 

The conditions (I)- (IV) are rather general, so that at any given order there are many 
functions satisfying these criteria. Concrete models of P^^^ (a) would typically contain some 
arbitrary parameters (this is the price we have to pay for getting rid of Landau singularity). 
One may restrict this freedom by fixing these parameters with help of some information 
coming from outside of perturbation theory. However, if in every order of perturbation 
expansion there would appear new, uncorrelated free parameters, the predictive power of 
such an approach would be very limited. For this reason we shall further restrict the class 
of possible functions (3^^\a). We shall require, that: 

(V) /3^^\a) should not contain free parameters specific to some particular order of the 
expansion (i.e. A^). 

In other words, the free parameters appearing in the modified /3-function should characterize 
the whole sequence of functions /3^^^(a), not just a model /3-function at some particular 
order. 



3.3 Concrete model of the modified couplant 



Nonpolynomial /^-functions that do not involve exponentially small contributions and which 
lead to the expressions for a(Q^) that are nonsingular for real positive have been con- 
sidered before ||l2l|, |l2|, |l2|, |l2|, |l2|, |l26|, |156|1. The simplest idea is to use the Fade 
approximant. In the NL order we then obtain 



1 



1 — ClO 



(3.7) 



Unfortunately, this expression becomes singular for real positive a, so there would not be 
much of an improvement in the behaviour of a. In higher orders the situation becomes 
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even worse — not only the Pade approximants develop singularities for real positive a, 
but also the locations of these singularities are scheme dependent. For example, the [1/1] 
approximant in the NNL order has the form: 

P^'Ha) = -U \ • (3.8) 

This shows that in order to obtain a satisfactory model /^-function one has to go beyond 
the straightforward Pade approximants or their simple modifications. 

Let us therefore construct a completely new model /3-function. Looking for a system- 
atic method to produce nonpolynomial functions with the correct weak-coupling expansion 
coefficients and appropriate asymptotic behavior we choose an approach inspired by the so 



called method of conformal mapping [|157| , [158 , 15G, 160, 161]. This method is essentially a 



convenient procedure for performing an analytic continuation of the series expansion beyond 
its radius of convergence and it has been a popular tool in resumming the divergent series 



via the Borel transform |159, 16C, |l6l|. (However, our use of the conformal mapping would 
be closer to the attempts to resum the divergent series by an order dependent mapping 
||75|] .) The basic idea of the mapping method is very simple: we take a function u(a), which 
maps the complex-a plane onto some neighborhood of u = in the complex-u plane and 
which has the following properties: (a) it is analytic in some neighborhood of a = 0, with 
the expansion around a = of the form u(a) = a + 0(a^); (b) it is positive and nonsingular 
for real positive o, with an inverse a{u); (c) in the limit a — > oo it goes to a real positive 
constant. In the following we shall use a very simple mapping 

uia) = (3.9) 

where ry is a real positive parameter, but of course other choices are possible. 

Let us now apply the mapping method to obtain j3{a) with the asymptotics correspond- 
ing to A; = 1 in the condition (III). We take the expression 

+ cia^ + C2a^ + ... + CAra^+^ 



and rewrite it as 



Ka — Ka + 0? + cia^ + C2a^ + ... + CAra^"'"^. 



Then we omit the first term 

—na + + cio? + C2a^ + ... + cncl^^'^ , 

and we substitute everywhere in this expression a = a(u); expanding the resulting function 
in powers of u, we obtain 

— KU + CqU^ + ClU^ + C2U^ + ... + CnU^'^'^ + ■■■ 

It is easy to see that the function 

/3W(a) = [na - Ku{a) + cou{af + ciu{af + C2u(a)^ + ... + CNu{a)'^+^] . (3.10) 
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does indeed satisfy the conditions (I)-(III) and (V) described in the previous subsection. 
For the concrete mapping (3^) we have 



Cq = 1 — rjK, ci = Ci+2ri — rfn, (3-11) 

and 

C2 = C2 + Sciry + 3ry2 - ry^K. (3.12) 

More explicitly, the modified /3-function that we propose to use has in the NNL order the 
form 



Ka , , a? 



Ko-— + {l-r]Hi)— (3.13) 

1 + rja (1 + r/o)^ 



{ci + 2r] - r] k)— -tt + (c2 + Scif] + 3ry - ry k 



(1 + ??a)3 ■ V- — / . " , ' (1 + r,aY 

It remains to verify that the function (|3.13| ) satisfies also the condition (IV). Expanding 
P^^^ (a) in terms of a we find 

/3W(a) -/3W(a) = ^ r] [Sci + 3rj - rj'^ k] a'^ + 0{a^), (3.14) 

and 

/3^^\a) - p^^\a) = ^r/ [4c2 + 6ci?7 + V - 7^3^] + 0{a^), (3.15) 

which shows that our simple model /3-function does indeed satisfy also the condition (IV), 
at least in the NL and NNL order. (This is to some extent a lucky coincidence — there 
does not seem to be any simple constructive method to satisfy the condition (IV).) 

Our candidate for the modified /3-function contains two free parameters: the parameter 
At, which determines the exponent in the low-Q^ behaviour of the couplant a(Q^), and the 
parameter rj, the inverse of which characterizes the range of values of a, for which the non- 
polynomial character of the (3^^^ (a) becomes essential. For any values of these parameters 
the function f3^^\a) is a viable replacement for P^^\a). However, by a proper choice of 
these parameters we may improve the quality of the modified perturbation expansion in 
low orders. As a first try, we shall assume in further calculations reported in the following 
subsections the value k = 2/6, which ensures ^ = 2 in the Equation ( |3.6D , implying l/Q'^ 
behaviour of a{Q'^) at low Q^, as is suggested by some theoretical approaches ||162|| . In a 
more precise approach one could consider adjusting this parameter in some range, together 
with the parameter rj. The procedure for choosing a preferred value for rj is described 
further in the text. 

The plot of this model /3 in NL and NNL order for a particular value of t] is shown in 
Figure 0. It is easy to guess from this figure, what the plots for smaller values of rj would 
look like, since in the limit r/ — > the function ^^^\a) coincides with the conventional 
/3-function. 

In order to determine the Q^-dependence of the modified couplant we have to calcu- 
late the function F^^)(a, C2, ctv, •••) defined by Equation ( |2.5|) , with j3^^\a) replaced by 
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Figure 7: The modified /3-function, as given by Equation (p. 131 ) with k = 2/b and rj = 4.1, in the 
NL order (dash-dotted line) and the NNL order for three values of C2 (soHd Hnes), compared with 
the conventional /^-function (dashed lines). 



f3^^\a). The relevant integral in the NL order may be calculated in closed form. The 
function F^^\a,r], k) in the Equation ( |A.3| ) takes the form: 



where 



F«( 



a,r],K) 



-E2 In [1 + (ci + 3?7)a + r/^Ka^] + 
ci + 3r? 



+E'i arctan 



-E-x arctan 



\/—ci — 6ci?7 — 9r/2 + Arj^K 



E2 
E3 



1 + ClK 

2k 

— ci — 3r/ + cf K + 3ciryK + 6r/^«; — 2r]'^K 



3,^2 



K\J —c\ — 6ci?7 — 9ry2 + 4773^ 
Strictly speaking, this expression is valid for 



(3.16) 



-cf - 6cir? - 9r/^ + 4?7^k > 0. 
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Figure 8: The modified couplant a{Q^), corresponding to the /^-function ( 3.13 ), in the NL order 
(dash-dotted fine) and NNL order (solid lines) for three values of C2, as given by Equation (A.3) 



with appropriate functions F^^^ , for 



Ag^=0.35GeV, 



2/& and 77 = 4.1. 



The expression for F'^^) (a, rj, k) for other values of the parameters is obtained by analytic 
continuation. The integrals in the NNL order are more complicated, but still may be 
handled by Mathematica. 

Examples of the plots of the modified couplant a{Q^) in the NL and NNL order for 
a particular value of rj are shown in Figure |8|; these plots should be compared with the 
corresponding plots for the conventional couplant, shown in Figure ||. We clearly see that 
the modified couplant is much less scheme and order dependent than the conventional 
couplant. 



4. Modified perturbation expansion 

4.1 Perturbative predictions with a modified couplant 

In the previous section we constructed a modified couplant a(Q^) that is free from Landau 
singularity and is much less RS dependent than the conventional couplant. The modified 
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Figure 9: The conventional couplant a(Q^) in the NNL order, as given by Equation (A^) with 
appropriate functions F'^'^\ ri — and = 0.35 GeV, for three values of C2. 



couplant is determined in the implicit way by the Equation (|2y 

^In-^— = rj°^^^ -n +ciln^ + i + ci In a + F(^)(a, C2, cat, r/, k), (4.1) 

^MS ^ " 

which is obtained by integrating (with an appropriate boundary condition) the modified 
renormalization group equation ( |2.2| ) 

0'^ = ^™©. m> 

with a modified generator f3^^^ given by Equation ( |3.13| ). 

We shall now use this new couplant to construct modified perturbative expressions 
for some physical quantities. We will show that the predictions obtained in the modified 
expansion are much less RS dependent than the predictions in the conventional expansion. 

The simplest way to obtain an improved perturbative expression for a physical quantity 
6 is to replace the conventional couplant a{Q'^) in the usual RG improved expansion for 
(5^^) by the modified couplant: 

5W(Q2) = ^(g2) ^ ~(g2) ^ ~2(g2) ^ ^^^^ a"" {Q^)] . (4.3) 
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The coefficients rj in the improved approximant are constrained by the requirement of the 
consistency with the conventional expansion. In fact, they may be taken to be the same as 
in the conventional approximant. This follows from the fact that the function P^^^ satisfies 
the condition (I) from Subsection 3.2, which in turn implies, that the difference between 
^(Af)(Q^) and a(N){Q'^) is of the order a^"*"^, so expanding the improved approximant 6^^^ 
given by ( |4.3| ) up to and including the order a^~^^ we recover the conventional N-th order 
approximant 6^^^ . This means, that we may characterize the RS dependence of the modified 
perturbative results using the same scheme parameters as in the case of the conventional 
approximants. 

In the following we shall usually omit the tildas, when there is no doubt that we are 
dealing with the modified couplant and the modified expression for the physical quantity. 

4.2 PMS scheme in the modified expansion 

In our analysis of the RS dependence of conventional approximants in Section 2 the PMS 
scheme played an important role. The PMS scheme may be defined also for the modified 
approximants, as we describe below. 

To obtain the modified PMS approximant in the NL order we apply the PMS condition 
( B.2 ) to the modified NL order approximant 6^^\ In this way we obtain an equation, which 



has the same general form as the corresponding equation for the the conventional expansion, 
except that the /?(^) is replaced by (3^^^: 



Solving this for ri we find 



a2 + (l + 2fia)-/3(i)(a) = (4.4) 



6a2 + 2/3(i)(a) 

n(a) = ,.5(r,' ■ 

4a;3(i)(a) 



For defined by Equation ( |3l3| ), with the mapping ( |3.9| ) and the coefficients ( p.ll| ), this 
implies 

^ -C1 + 77^(3 -7?K)a + r?3a2 
' 2 [1 + (ci + 3r?)a + ??3Ka2] • ^ > 

Inserting this into the Equation (^) for the NL order modified couplant (with the function 
given by (|3.16D ) we obtain certain transcendental equation for a; solving this equation 

and inserting the relevant values of a and fi into the expression (p|) for (5(1) we then obtain 

the NL order PMS prediction for 6 in the modified expansion. 

For larger values of we may obtain an approximate expression for fi in the form of 

expansion in powers of a. It is interesting that the leading order term in this expansion is 

independent of the concrete form of /^^^^ and coincides with the value ( |B.4| ) obtained in the 

case of the conventional NL order PMS approximant: 

n = -| + 0(a), (4.7) 

as may be easily verified by inserting P^^\a) = P^^\a) + 0{a'^) into (|4.5D . This approximate 
solution does not have much practical value for our considerations, since we are interested 
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in the region of moderate . However, it is interesting from the theoretical point of view 
and may serve as a cross check for various numerical procedures. 

The PMS conditions (|2.7D and (|2.8D for the NNL order approximant in the modified 
expansion generate equations, which in the general appearance are similar to the equations 
obtained for the conventional expansion, although of course they are more complicated. 
The partial derivative needed for the PMS equation (2/7) has the form: 

Ap) = ^2 ^ + 2rx)c? + (1 + 2ria + ■iT2a^)\W\a), (4.8) 
or\ 

where we used the notation 6^"^^ to emphasize that we are dealing with the modified expres- 
sion for 6^"^^ A more explicit form of this derivative is not very instructive, except for the 
leading term in the expansion in powers of a, which is important for further considerations: 

Ap) = _(3p2 _ 2c2 + 5ciri + 3rl)a^ + 0{a^), (4.9) 

It is interesting that this term is independent of the concrete form of (3^"^^ and has the same 
form as the corresponding term in the case of the conventional expansion (Equations ( [B.7D 
and ( [B.8| )). This may be easily verified by inserting (3^'^\a) = P^'^\a) + O(a^) into the 
Equation (gj). 

The partial derivative needed for the PMS equation (|2.8| ) has the form 



dc2 Jo (/5(2 



(/3(2))2 dc2 ' 



(4.10) 



which coincides with the Equation ( [B.9D , except that /J^^) is replaced by Again, the 

explicit form of this derivative is not very illuminating, except for the leading order term in 
the expansion in powers of a. This term depends on the leading order difference between 
and which we shall write in the following form: 



/3(2)(a)=/?(2)_^d3a5 + 0(a6). (4.11) 



Using this notation we obtain: 

In the concrete case of given by ( |3.13| ), (p^), ( |3.11] ) and ( p.l2| ) we have 



dco 



(2ri -2r/)a* + 0(a^). (4.13) 



The PMS prediction for <5(2) in the modified expansion is obtained in the same way, as 
in the conventional expansion: solving the PMS equations (|B.8| ) and ([B.13| ) and the implicit 
equation for the NNL order couplant (|2.4D (with function F(2) appropriate for the modified 
expansion) for the chosen value of we obtain the parameters a, fi and C2 singled out 
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by the PMS method; inserting these values into the expression (^) for (5*^^^ we obtain the 
modified PMS prediction for 5^'^^ at this value of Q^. 

Similarly as in the case of the conventional PMS approximant, for small a (i.e. large 
(5^) we may obtain an approximate solutions of the PMS equations in the leading order 



in the expansion in powers of a, using the formulas ( [4 .91 ) and ( 4.12 ). Taking into account 



( 4.13 ) we immediately obtain from Equation ( p.SP that 

7=1= r] + 0{a). (4.14) 



Inserting this into the Equation (|2.7D and taking into account (|4.9| ) we then find: 

C2 = ^P2 + ^cir] + ^ri^ + 0{a). (4.15) 

In the limit ?7 — > these approximate solutions coincide with the corresponding approxi- 
mate expressions for the NNL order PMS parameters in the conventional expansion (see 
Appendix B). As we already mentioned in the case of the NL order approximant, these 
approximate solutions are of little practical significance at moderate Q^, where a is rather 
large, but they serve as a useful cross check for various numerical routines. They are also 
interesting from the theoretical point of view, because they show that — at least in the 
weak coupling limit — the solutions of the PMS equations for the modified expansion indeed 
exist and are unique. 

Let us note that the PMS scheme for the modified expansion stands theoretically on 
a better footing than in the case of the conventional expansion, because in the modified 
expansion the perturbative approximants remain finite for all scheme parameters and all 
positive values of Q^. 

4.3 Choosing the value of r/ 

By introducing the modified couplant we obtained an improved perturbation expansion 
that does not suffer from the Landau singularity and (as we shall see in the following 
subsections) is much less sensitive to the choice of the renormalization scheme. However, 
in this new expansion the finite order perturbative predictions acquire a dependence on 
the parameter ry. (They depend also on the parameter k, but to simplify the analysis 
we have initially assumed k = 2/6.) Erom the way we introduced the modified couplant 
it is clear that the sum of the perturbation series should be independent of r/, but this 
statement requires some comments. Eirst of all, the perturbation expansion in QCD is 
divergent, so in order to achieve the r/-independence we would have to use an appropriate 
summation method. As mentioned in Subsection 2.2, there is a possibility, that in the 
case of QCD such a summation method may be provided by the PMS approach, which 
motivates us to give some preference to the modified PMS approximants when considering 
phenomenological applications of the the modified expansion. Secondly, in fact we do 
expect some ry-dependence even if the perturbation series is resummed to all orders. This 
is because for ?7 = 0, i.e. in the conventional expansion, there are whole classes of schemes, 
for which resummation of the perturbation series for the physical quantity cannot give a 
satisfactory result, because it is unlikely that it would remove the Landau singularity in 
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the Q^-dependence of the couplant; one example of such a scheme is the so called 't Hooft 
scheme, in which all the /3-function coefficients beyond the NL order are chosen to be 
identically equal to zero: Cj = for i > 2. On the other hand, if r/ / 0, then perturbative 
predictions are finite for all positive values of Q^, and one may reasonably expect that the 
resummed series would give an expression that does not suffer from the Landau singularity 
problem. 

Even if we use some resummation procedure, the low order approximants do depend 
on r/. It should be stressed that such dependence may be in fact seen as an advantage of 
the modified expansion, because it may be used to improve the accuracy of the predictions 
by incorporating some information outside of perturbation theory. This may be achieved 
for example by adjusting i] to match some phenomenological or nonperturbative results for 
some important physical quantity. Ideally one would want to use some results obtained from 
first principles, for example in the lattice approach, but this is a complicated subject, so as 
a first try we use in this article a much simpler condition that relies on a phenomenological 
formula for 6y |12E, 130 1 that has had some success in correlating the experimental data 
for heavy quarkonia. This formula has the form of an implicit equation for dy (denoted as 
(^bgt): 



^ 1 MS,V , ^ 

— hi 77 — = r, H — 

2 |'A^)2 1 2 
^ MS ^ 



( b6BGT ) 



+ Cl 



. 26 ,A(5bgt^ 
ln--7E-Ei(^^, 



(4.16) 



where 7e is the Euler constant, Ei(x) denotes the exponential integral function, b and ci 
are the usual renormalization group coefficients for Uf = 3 and rj^^'^ is the first expansion 
coefficient for 6y in the MS scheme. The recommended values of the parameters in this 
expression are A = 24 and = 500 MeV. The expression ( [1.16|) is very convenient 

from our point of view, because it refers directly to the space and by construction it 
is consistent with the NL order perturbative asymptotic behaviour of (5v((5^) at large Q^. 
We propose to choose rj in such a way that at moderate the modified predictions for 
6y would match the phenomenological expression ( 4.16 ) as closely as possible. (It should 
be noted that a generalization of the expression ( [4.16| ) was recently discussed in [132|. 
However, authors of |132| incorporate in their formula a large NNL order correction and 
use as a constraint in the fit a rather large value of as(M|), so we decided to use the older 
expression (P^) . In any case, we have verified that the fitted value of r] is not affected in 
a significant way.) 

The concrete procedure that we use to fix r/ is the following: we assume that the mod- 
ified predictions for coincide with the value given by the phenomenological expression 
( |4l6|) at Q2 = 9 GeV^, i.e. at the upper boundary of the n/ = 3 region (which is achieved by 
adjusting the parameter Ay^ in the expression for 5y) and we adjust rj so as to obtain the 
best possible agreement between these two expressions at lower Q^, down to = 1 GeV^. 
To perform this matching procedure we use the predictions for dy in the PMS scheme. (We 
justified our preference for this scheme above.) Putting aside the problems of resummation, 
we could of course fix r] using predictions in other renormalization schemes, and the result- 
ing value would of course depend to some extent on the scheme; however, this dependence 
would be very small, because the scheme dependence of the modified predictions is small 
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(as we shall see), and the difference to large extent would cancel away, if we would calculate 
the predictions for other physical quantities in the same renormalization scheme. 

The curves illustrating the relative deviation of the modified PMS predictions for 
in the NNL order from the phenomenological expression ( [4.16| ) are shown in Figure [l^ for 
three values of t]. It appears that the best fit is obtained for 77 « 4.1 — for this value of rj 
the relative deviation is less than 1% down to = 1 GeV^. This is the value that we shall 
use in further analysis. 
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Figure 10: The relative difference between the modified NNL order PMS prediction for Sy and the 
value obtained from the phenomenological expression (4.16), as a function of Q^, for three values 
of 77. In this calculation the parameter in the expression for Sy is adjusted in such a way that 
it coincides with (5bgt at = 9 GeV^ . 



Fixing the value of rj we should take into account the fact that the phenomenological 
expression (|1|) has hmited precision. In order to estimate, how this might affect the 
preferred value of rj we repeated the fitting procedure, using = 400 MeV in the phe- 
nomenological expression (|1|). (This value has been quoted in |129| ] as a lower limit on 
this effective parameter.) We found that this change has only slight effect on the fitted 
value of ry — the general picture is unchanged and a good match is obtained for rj = 4.0. 
This indicates that the uncertainty in the fitted value of r/ arising from the uncertainty in 
the phenomenological expression is of the order of ±0.1. 
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4.4 6y and (5gls in the modified expansion 

Having fixed r] we may now verify, whetlier tlie use of the modified couplant indeed results 
in a better stability of tlie predictions. In Figure |l^ we show the modified predictions for 
6y as a function of ri, for = 3GeV^ and rj = 4.1. In Figure |l^ we show the modified 
predictions for the same quantity, in the NNL approximation, as a function of Q^, for several 



values of ri and C2. In Figure and Figure |lj we show the corresponding plots for (5gls 



Comparing these figures with Figures HI § and respectively, we may clearly see that 
the scheme dependence of the modified predictions is substantially smaller (at least for this 
value of r]) than the scheme dependence of the predictions obtained from the conventional 
expansion. This means that our initial hypothesis concerning the origin of the strong RS 
dependence of conventional approximants (that it is due to strong RS dependence of the 
conventional couplant itself) was correct. 

Let us note, that strictly speaking, in our comparison of the RS dependence of the 
modified and conventional perturbative predictions the modified predictions should be 
plotted with a slightly different value of Aj^. Indeed, if we take the world average of 



as{M'^) = 0.1182 ± 0.0027 found in [|T0|1 and calculate the value of A^^, for which the 



same value is obtained for the modified strong coupling constant in the MS scheme at the 
Mz scale (calculated with rj = 4.1, using the matching procedure described in detail in 
Section 6), we get A^ = 0.37 GeV±Q;Q4. However, changing from A^ = 0.35 GeV to 
0.37 GeV does not affect the plots shown in this subsection in any significant way, so for 

(3) 

the sake of simplicity we use A^ = 0.35 GeV in all the plots throughout this article. 
4.5 PMS predictions in the modified expansion 

In Subsections 2.2 and 4.3 we mentioned some arguments in favor of the PMS scheme, so it 
is of some interest to discuss the PMS predictions for 6y and (5gls in more detail. Looking 
at the Figures 4.2 and 4.3, which show the reduced RS dependence of the predictions in 
the modified expansion, we make an interesting observation concerning the predictions in 
the PMS scheme. In the case of the conventional expansion at moderate the variation 
of predictions for finite (i.e. not infinitesimal) deviations of the scheme parameters from 
the PMS values is quite strong, as may be seen in Figures 2.1 and 2.5, which could raise 
doubts about validity of this prescription to choose the scheme. In the case of the modified 
expansion the scheme dependence in the vicinity of PMS parameters much weaker, simply 
because the modified predictions are generally much less sensitive to the choice of the 
scheme parameters. This further supports our conclusion at the end of Subsection 4.2 
that the PMS scheme in the modified expansion stands on a better footing than in the 
conventional expansion. 



In Figure 15 we show the curves representing the modified NL and NNL order predic- 
tions for 6y in the PMS scheme, compared with the conventional PMS predictions for this 
quantity. Strictly speaking, the NNL order curve for 5y is not a pure prediction, because 
we used this quantity to fix the value of the parameter rj (although this fitting had a rather 
indirect character — we fitted only the energy dependence, not the numerical values of the 
predictions themselves). However, having done that, we are in position to make some gen- 
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Figure 11: Sy &t Q"^ = 3 GeV^ (for n/ = 3) as a function of ri , for several values of C2 , obtained 



in the NNL order in the modified perturbation expansion with the /^-function (3.13) and i] — 4.1. 
Dashed line indicates the NL order prediction. 



uine improved predictions for any other physical quantity, for example Sqls- In Figure 16 
we show the NL and NNL order modified predictions for Jgls in the PMS scheme, compared 
with the PMS predictions for this quantity obtained in the conventional expansion. 

One characteristic feature shown by these figures is that the PMS predictions in the 
modified expansion lie below the PMS predictions in the conventional expansion, both in the 
NL and NNL order. This effect is perhaps not very surprising, given the way we constructed 
the modified couplant, but a moment of thought shows, that it is also nontrivial, since we 
are dealing with "floating" optimal scheme parameters that in the modified expansion have 
different value than in conventional expansion, as shown for example by the Equations ( 4.14 ) 
and (4.15). We discuss this issue in greater detail in Subsection 6.1. 



5. Modified expansion with a different couplant 
5.1 Alternative definitions of the modified couplant 

Before we apply the modified expansion to any phenomenological problems, we have to 
address an important question, to what extent the predictions obtained in the modified 
expansion depend on the exact form of the chosen modified /3-function, since the set of 
functions satisfying the conditions (I)-(V) from Subsection 3.2 is still rather large. In order 
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Figure 12: (Sy as a function of (for Uf = 3), obtained in the NNL order in the modified 
perturbation expansion with the /3-function ( 3.13 ) and 77 — 4.1 in several renormaHzation schemes, 
including the MS scheme (dash-dotted line) and the PMS scheme (soHd line). 



to investigate this problem in a quantitative way let us consider an alternative class of 
the modified /3-functions, resulting from a different application of the mapping method, as 
described in Subsection 3.3. We start with the expression 

a^-''{l + cia + ... + CNa^), (5.1) 

where A; < 1; we then substitute a{u), extract the factor u^~^ in front of the expression 
(because for general k it would give an essential singularity) and expand the rest in the 
powers of u. This gives 

u^^^{l + ciu + ... + CNU^ + .:) (5.2) 

where Cn are appropriately modified expansion coefficients. As a modified /3-function we 
now take the following expression: 

/3W(a) = -^-a^ {u{a)f-'' [l + cin(a) + ... + c^(n(a))^] (5.3) 

To make further progress we have to specify explicitly the form of u{a). If we use the 
same mapping as previously (Equation (^)), we find 

ci = ci + (2 - fc)ry, 02 = 02 + ^r/(3 - k) [(2 - k)r] + 2ci] . (5.4) 
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Figure 13: (5gls at Q2 = 3GeV^ (for Uf = 3), as a function of ri , for several values of C2 , obtained 
in the NNL order in the modified perturbation expansion with /3-function (3.13) and rj = 4.1. 
Dashed line indicates the NL order prediction. 



In this way for any chosen /c < 1 we obtain a very simple modification of P^^\a) that 



in contrast to (|3.13|) contains only one free parameter r/. It is obvious from the above 
construction that /3-functions constructed in this way indeed satisfy the conditions (I)-(III) 
and (V), formulated in Subsection 3.2. 

Re-expanding /3(a) in powers of a we find in the NL order 

/3(i)(a) _ = ^ [2ci(3 - k)7] + (3 - k){2 - k)r,^] + 0{a^), (5.5) 

and in the NNL order 

/3(2)(a)-/3(2)(a) = ^{4.-k)r] [6c2 + 3ci{3 - k)r] + {3 - k){2 - k)r]^] a^+0{a^) (5.6) 

As we see, the condition (IV) is also satisfied, at least up to the NNL order. 

There is a price we must pay for the simplicity of these models for P^^'^ (a), namely that 
the coefficient in the large-a behaviour of (3^^^ cannot be freely adjusted and is uniquely 
fixed by rj and the coefficients of the small-a expansion. This is however perfectly consistent 
with our general approach, since our goal is only to improve the predictions at moderate 
Q^, and from this point of view the exact behaviour of the couplant at very low is not 
very important. 
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Figure 14: (5gls as a function of (for n/ = 3), obtained in the NNL order in the modified pertur- 
bation expansion with /^-function ( 3.13 ) and 77 = 4.1, in several renormalization schemes, including 
the PMS scheme (solid line). The curve corresponding to the MS scheme is indistinguishable from 
the PMS curve in the scale of this figure. 



For concrete numerical calculations we choose j3^^\a) of the form ( ^I^ ) with A; = 
(the /3-function of this form has been first discussed by the present author in |E3|): 



a = - 



2(l + r?a)2 



1 + (ci + 2r]) 



1 -\- r]a 



+ (c2 + 3ci7y + 37]'^ 



1 + 7]a 



(5.7) 



The plot of this /3-function in the NL and NNL order for a particular value of t] is shown 
in Figure |l^ (in this plot we assumed rj = 3.8, which is justified in Subsection 5.2). It is 
interesting that although this /3-function has different large-a asymptotics than the function 
given by Equation ( ^.13 ), the qualitative behaviour of these functions in the range of a which 
is important at moderate is quite similar. 

The corresponding function F^^^ in the Equation ( ^.4] ), determining the Q^-dependence 
of the couplant, has in this case in the NL order a very simple form: 



F«(a,ci,r/) 



(ci + 2rjY 
(ci + SrjY^ 



In [1 + (ci + 3r])a] 



V 

ci + 3r] 



(5.5 
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Figure 15: The NL and NNL order PMS predictions for Sv, as a function of Q^, obtained in the 
modified expansion with /3-function ( 3.13 ) and 77 = 4.1 (sohd fines), and in the conventional RG 
improved expansion (dash-dotted fines). 



The expression for F^'^^ in the NNL order for this /3- function is lengthy, but more manageable 
than the corresponding formula for the function ( ^.13|) : 

-E2ln[l + a (ci+477)+a^ (c2 + 4ci r/ + Gr/^)] + 

ci + 2 77 \ 



where 



+ E3 arctan 



y'-ci+4c2 + 8cii] + 8r]^ 



/ ci + 2?? + a (2c2 + 7cir? + 8r/2) \ 

-E/3 arctan , (5.9) 

Ul + ar]) V-c^ + 4c2 + 8ciry + 8?72y 



^ _ ci + 8 cf C2 T? + (16 cf + 12 ci C2) r?^ 
^~ 2 (C2 +4cir? + 6r?2)2 

(48 cf + 4 C2) + 51 ci r/'' + 20 

2 (c2 + 4ci7? + 6 7?2)2 
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Figure 16: The NL and NNL order PMS predictions for (5gls, as a function of Q^, obtained in 



the modified expansion with /?-function ( 3.13 ) and 77 = 4.1 (soUd lines), and in the conventional 



(5.11) 



RG improved expansion (dash-dotted lines), 
and 

^ _ (cf-2c2) cl + Aci (2cf-5c2) C2 r? + 4 (4 cf - 13 C2 - 6 cj) 7?^ 
^ (c2 + 4ci77 + 6r/2)^ ^_c^ + 4c2 + 8cir/ + 8r/2 

4 (4cf + 37C1C2) 77^ + 9 (I9cf + IOC2) + 232 ci + 92 
(c2 + 4cir? + 6r?2)2 y^-c^ + 4 C2 + 8 ci r? + 8 7?2 

Strictly speaking, this form is valid for a, C2 and satisfying the condition 

-c? + 4c2 + 8cir/ + 87?2 >0, 

which also guarantees that 

1 + (ci + 4r?)a + (c2 + 4cir/ + 6??^)a^ > 0. 

For other values of the parameters C2 and rj the function F^^) is obtained via analytic 
continuation. 

Similarly as in the previous case we define the modified expansion for a physical quantity 
6 by replacing the conventional couplant in the series expansion for 6 with the modified 
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Figure 17: Alternative model of a modified /3-function, as given by Equation (|5.7D with ry = 3.8, 
in the NL order (dash-dotted line) and the NNL order for three values of C2 (solid lines) , compared 
with the conventional /3-function (dashed lines). 



couplant determined by (|5.7| ). The PMS approximants are also defined in a similar way as 
in Subsection 4.2. The NL order approximant is defined by the equation (|4.4| ), where (3^^^^ 



(5.12) 



is now given by ( |5.7|) . Solving this equation with respect to ri we obtain: 

" 2(1 + (ci + 3T])a) 

The PMS prediction for 5^^^ in this modified expansion is then obtained in the same way 
as in Subsection 4.2. 

The NNL order PMS approximant for this model of the modified couplant is defined 

and 



by the Equations (2/7) and (|2.8|), with the relevant partial derivatives given by 
dilOD , where /J^^) is now given by (|5.7|). Similarly as in the previous case, these equations 
are too cumbersome to be reproduced here in the full form, but it is easy to obtain an 



explicit approximate solution for them in the leading order of expansion in a. Using ( 4.12 ) 
and taking into account the form of the coefficient d^, which follows from the expansion 
dp), we find: 



dc2 



(2ri -2r/)a^ + 0(a^). 



(5.13) 



Somewhat surprisingly, this coincides with the expression (fl.l3|) , obtained for the previous 
model of the modified couplant. This means that in the leading order in a the approximate 
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values of the PMS parameters for the coupling defined by (5/7) are the same as the ap 



proximate PMS parameters ([4.14|) and ( 4.15 ), obtained for the couplant defined by ( 3.13 ). 



However, looking at the Equation (|4.12 ) for the approximate value of the derivative of 5^^^ 



over C2 and taking into account the form of the relevant coefficient in the Equation (5.6) 
we immediately see that this is just a coincidence, i.e. for other values of k in the modified 
/3- function ( |5.3D the approximate values of PMS parameters would be different. As men- 
tioned previously, this approximate solution does not have much practical importance for 
our discussion, but it shows that also for this type of the modified couplant the solutions 
of the PMS equations in the NNL order exist and are unique, at least in the weak coupling 
limit. 

5.2 Choosing t] for the alternative modified couplant 

Similarly as in the previous case, the modified predictions in low orders of the expansion 
depend on the value of the parameter r/, and in order to obtain meaningful predictions we 
have to fix somehow this parameter. We use exactly the same procedure as described in 
Subsection 4.3. The value giving the best fit appears to be approximately ry = 3.8; for this 
value the relative deviation is within 1% down to = 1 GeV^. 

5.3 6y and Sgls for alternative modified couplant 

In Figure 18 we show at = 3 GeV^, as a function of ri, for several values of C2, obtained 



with an alternative couplant defined by the /3-function (|5.7|) . As we see, this figure is very 
similar to the corresponding Figure |ll], obtained with the previous model of the /3-function. 
The same is true with respect to ri dependence of the NNL order predictions for (^gls 
at fixed and the dependence of the NNL order predictions for both quantities for 
various (fixed) scheme parameters. 



In Figure 19 we show the PMS predictions for 5gls, obtained in the modified expansion 



with the couplant defined by Equation {f>-7\j. This figure should be compared with Figure [16 
representing PMS predictions for the same quantity, obtained with the modified couplant 
defined by the Equation ( 3.13| ). We find that the NNL order predictions obtained in both 



cases are practically identical. However, the NL order modified predictions do not agree so 
well, which is not surprising; they are nevertheless consistent, provided we take into account 
the fact, that one should assign to those predictions the theoretical error of the order of 
the difference between NL and NNL order predictions in each case. Let us note that the 
couplant defined by Equation ( ^.71 ) leads to larger NL/NNL difference, which means that 
the modified expansion in this case is less precise. 

The consistency of the predictions obtained with different modified couplants is rather 
encouraging, because it shows, that despite some freedom in choosing the modified /3- 
function our approach leads to consistent physical predictions even in low orders of per- 
turbation expansion. We see two factors that might be responsible for this "universality" 
of the predictions. One is that we have an adjustable parameter rj, which to some extent 
absorbs possible differences. Second factor is the use of the PMS approximant for physical 
predictions: using PMS we no longer choose the RS parameters by hand — the values of 
these parameters "adjust themselves" in an "automatic" way; presumably the exact behavior 
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r igure 18: Sy at Q2 = 3 QeV^ (for n/ = 3), as a function of ri , for several values of C2 , obtained 
using the alternative modified couplant defined by Equation (5.7) with 77 = 3.8. Dashed Hue 
indicates the NL order prediction. 



of the running couplant at low energies is not so important, as long as it stays finite for all 
positive Q^. This gives additional motivation for using PMS in the modified approach. 

We might expect that this behaviour would be a general feature of the modified per- 
turbation expansion, i.e. concretely if we use the NNL order expression for some physical 
quantity to fix the adjustable parameters in /?, then the NNL order predictions for all other 
physical quantities should not be significantly altered by changing the form of the modified 
/^-function, at least when the predictions are evaluated in the PMS scheme. However, the 
differences between predictions in successive orders for the same physical quantity may vary, 
depending on /3. 



6. Consequences for phenomenology 

6.1 Reduced Q^-dependence of the modified predictions 

Our analysis of the perturbative predictions for 5y and 5gls has shown an interesting 
effect: the PMS predictions in the modified expansion lie below the PMS predictions in the 
conventional expansion. This effect is of course not unexpected, since the initial motivation 
for our construction of the modified expansion was exactly to reduce the rate of growth of 
the couplant with decreasing in some schemes and to eliminate the singularity at positive 
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Figure 19: The NL and NNL order PMS predictions for (5gls, as a function of Q^, obtained with 
the alternative modified couplant defined by Equation (5/7) with rj = 3.8 (solid lines), compared 
with the conventional PMS predictions (dash-dotted lines). 



Q^. On the other hand, this effect is not entirely trivial, because there is no immediate 
relation between the low-Q^ behaviour of the couplant in any particular scheme and the 
low-Q^ behaviour of the PMS approximants: firstly, in some renormalization schemes the 
modified predictions lie above the conventional predictions, as may be seen for example in 
Figure |ll|, and secondly, the PMS scheme parameters for the modified approximants are 
different from the PMS parameters of the conventional approximants. If such a pattern 
would be observed for all physical quantities, it might have important consequences for 
phenomenology. It is therefore of some interest to investigate, whether we may support our 
observation by some general arguments. 

An ideal solution would be to give some rigorous bounds on the PMS predictions in 
the modified expansion, but this seems to be a rather difficult problem. We shall therefore 
adopt a less ambitious approach: in order to compare the modified PMS predictions with 
the conventional PMS predictions in a given order we shall expand both expressions in 
terms of a new couplant ao, defined by the implicit equation: 

b , Q'^ (0)MS ,1^,1,1 2 fa i\ 

-ln-2— =r}' +ciln-H h ci In oq - CiOq, (6.1) 
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where rj'^^'^^ denotes as usual the first expansion coefficient for the considered physical 
quantity in the MS scheme. The usefulness of the couplant oq lies in the fact that it is 
nonsingular for all real positive and at the same time it has correct weak coupling 
behaviour up to NL order. It satisfies the renormalization group equation of the form 



2 dao 

2 1 — Ci Oq + cf Oq 



= ^^^^ (6-2) 



For the A^-th order modified PMS approximants the first correction terms in the ex- 
pansion in terms of ao must coincide with the corresponding terms in the expansion of the 
conventional PMS approximants (this follows from the general properties of our modified 
expansion), so in order to make an estimate of the difference between these two types of 
approximants we have to compare the first correction terms beyond the formal order of the 
approximants. 

Let us first consider the modified NL order PMS approximant, defined by the Equation 
(|4.6| ). We want to calculate the ©(oq) term in the expansion of this approximant in powers 



of Oq. To this end we insert (|4.6|) into (2.4) and we assume that the modified couplant 



a((5^) in the PMS scheme is related to ao(Q^) via the relation 

o = ao [l + ^1 ao + Ai ag + A3 aj] + ...] . (6.3) 
Making use of the defining equation ( |6.lD for ao((5^) we find: 



^1 = 2^1^ ^2 = |cf--cir?--77^ + -K?7^ (6.4) 

Taking this into account and expanding the expression for ^^^^ in terms of ao we obtain: 

5(1) = ao + Q - 3ci - 3r/2 + Kjf^ a|] + O(a^). (6.5) 

The relevant expansion for the conventional NL order PMS approximant is obtained from 
this result by setting r/ = 0. We see therefore that the difference between the NL order 
PMS approximants in the modified expansion and the conventional expansion is given by 

41 - ^^o). = (-3ci r? - 37?2 + nrf) aj] + O(a^). (6.6) 

In our calculations we assumed k = 2/6 = 4/9; for this value of k the coefficient of the 
leading term in this difference is negative for rj < 8.21. This shows, that the modified NL 
order PMS predictions would indeed lie below the conventional PMS predictions for all 
physical quantities that may be written in the form (|A.l ), at least for small ao. 



In a similar way we obtain the expansion of the NNL order modified PMS predictions 
up to the ©(cq) term. We first look for the solutions of the system of two equations 

d 

— 5(^)(a((3^,ri,C2),ri,C2) |^,=f,,c2=c2= 0, (6.7) 

and 

d 

— 5(^)(a(Q^,ri,C2),ri,C2) |ri=fi,c2=c2= 0, (6.8) 
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in the form of the power series in oq. We find: 



fi = 7] + Biao + B2al + ..., (6.9) 



where 



and 



where 



o 1 3 5 2 

27 29 3 

B2 = -P2r]-—ci7j^ -—r]^ + -Kr]\ (6.10) 



3 5 3 
C2 = ^P2 + ^cir] + -1]'^ + Diao + :., (6.11) 

^1 = ^ ci p2 + ^ c? r? + ^ ci r?2 _ ^3 ^ 1 ^^4_ ;^2) 



We then insert these values of fi and C2 into the defining equation ( |4.1| ) for the NNL order 
modified PMS couplant a, and we look for the solution for o in the form of the expansion 
( |6.3| ). Taking into account the defining equation for oq we now find: 

Ai = -ri 

5 2 

^2 = P2 + - r/ 

1 s 11 57 2 89 o 5 4 

^3 = 2Ci-yP2??-yCi7?-— r?+-Kr/ (6.13) 

Inserting the expressions for fi, C2 and a into the expression for 6^"^^ we find: 

5(2) = + p2 ag + Q c? - 3p2 r? - y ci r,^ - Sr?^ + i ^r/^^ + O(a^). (6.14) 

Again, the relevant expansion for the conventional NNL order PMS approximant is obtained 
from this expression by setting r/ = 0. We see therefore that the difference between the 
NNL order PMS predictions in the modified expansion and the conventional expansion is 
equal to 

€L - ^^l = (-3P2?? - y ci r?2 - 3r?3 + i Kr?^^ + O(a^). (6.15) 

For K = 2/6 the coefficient in the leading term is negative for rj < 16.40. This shows that 
(at least for small ao) also in the NNL order the modified PMS predictions would lie below 
the conventional PMS predictions for all physical quantities that may be expressed in the 
form ( A.l ). 

6.2 Some fits for the GLS sum rule 

The fact that the predictions for physical quantities obtained in the modified expansion 
lie below the predictions obtained in conventional expansion and evolve less rapidly with 
may have interesting consequences for phenomenology. To illustrate this, let us again 
consider the GLS sum rule and let us see, how the use of the modified expansion affects 
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the results of the fits to experimental data. We shall rely on the result reported in [107|, 
where it was found that measurements of deep inelastic neutrino-nucleon scattering for 
1< Q2 < 15 GeV^ imply 



af ^(3 GeV^) = 0.23 ± 0.035(stat) ± 0.05(sys). 



(6.16) 



Using the = 3 NNL order formula for (^gls we may interpret this result as 5^Lg(3 GeV^) = 
0.131 ±0.40, where we have added the statistical and systematic errors in quadrature. This 
would be the starting point of our phenomenological analysis. (To be precise, one should 
take into account at each value of the contributions from different values of nj, as 



discussed in |l 1 1]] , but this makes the whole discussion much more complicated. In order to 
see clearly possible effects associated with the use of the modified couplant instead of the 

(3) 

conventional couplant we limit our discussion to the case of n/ = 3.) We fit to this 
number using various theoretical expressions for (5gls, and then we convert the result into 
the corresponding value of a^^(M|). To relate the values of the couplant a = OLsji^ for 



different numbers of quark flavors we use the matching formula |162 , 164, |165| , 166, |167| , [168| 



«(A^^^/) = a(Ai^n/ - 1) 



l + iln-^a(^2^n.-l) + 
3 vriq 



+ : 



In^ii + ^ln^ 

mq 4 ruq 



11 



a^{n^,nf - 1) 



(6.17) 



where niq is the running mass of the decoupled quark, evaluated at the scale = niq. (For 
evolution of the couplant in the NL order we use the above formula restricted to the first 
two terms.) We use nic = 1.3 GeV and mj, = 4.3 GeV and choose as the matching points 
IJ? = Aniq. We use the same matching formula both in the conventional and the modified 
expansion (which is justified, because the difference between the modified perturbative 
expression and the conventional perturbative expression at given order is formally of higher 
order) . 

As a consistency check of our approach we first perform the fit in the MS scheme. 
Using the conventional NNL order expression we find as(Mf) = 0.1126ibQ;Q5i6- (We keep 
an extra digit in the quoted numbers and present asymmetric errors in order to minimize 
the roundoff errors and illustrate the sensitivity of the fitted values in a better way.) This 



number is slightly smaller than the result quoted in the original paper |107|; the difference 
may be presumably attributed to the different method of matching the couplant across 
the thresholds and/or different way of calculating the Q^-evolution of the couplant. If 
we use the conventional NL order approximant, we find as{M'^) 



0.1162±[j:[]?|^, which 



in turn is slightly higher than the value quoted in ||107|| . Note that both numbers are 
below the world average as (M|) = 0.1176 ±0.0020 obtained in and the world average 



as(M^) = 0.1182±0.0027 calculated in |105(| , and in the case of the NNL order approximant 
the difference is significant. 

Now if we perform the same fit in the MS scheme, using the modified perturbative 
expression with the coupling defined by the Equation ( ^.13 ) with = 4.1, we find in the 
NNL order as(M|) = 0.1147±[];[J?^^. (This is the modified coupling at the scale fi"^ 



Ml 
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and we used the modified renormalization group equation to evolve from /x^ = 3 GeV^ 
to //^ = M|; it is meaningful to compare this value with as{M^) from other sources, 
because it is exactly this value that we would insert for example into the expression for 
the QCD corrections to T{Z hadrons), if we would want to make a direct comparison 
with experimental data at fi^ = M^.) When the NL order approximant is used, the result 
is as(M|) = 0.1224ibQ;Q}54. As we see, if the modified expansion is used in the fit, then 
somewhat higher values of as(M|) are obtained. 

If we perform the same fits using the conventional PMS approximants, we obtain 
a,(M|) = 0.1097±[]:[J?g| in the NNL order and as(M|) = 0.1099 ±[]:[J?g| in the NL or- 
der. Note that these values are smaller than values obtained from the fit in the MS scheme 
and much smaller than the world average quoted above. However, if we use the PMS ap- 
proximants in the modified expansion with the /3-function ( ^.13 ) and rj = 4.1, we obtain 
a,(M|) = 0.1150±[]:g?|^ in the NNL order and ds(M|) = 0.1156±[}:[];| in the NL order. 
Similarly as in the case of the MS scheme, the use of the modified expansion gives higher 
values of as(M^), that is close to the present world average. In the case of the PMS ap- 
proximants this effect is much more pronounced than in the case of the approximants in 
the MS scheme. Interestingly, the difference between the NL and NNL order results, which 
may be taken as an estimate of the theoretical error, is in this case much smaller in the 
PMS scheme than in the MS scheme, both in the conventional and the modified expansion. 

Both in the MS scheme and the PMS scheme we find that the value of Q!s(M|) arising 
from the fit to the experimental data is higher when the modified expansion is used. This 
is a welcome effect, because it makes the low and high energy determinations of the strong 
coupling parameter more consistent. (The world averages for as{M'^) quoted above are 
dominated by measurements at higher energies, corresponding to Uf = 4,5.) The shift in 
the central value is not very big in comparison with the experimental error for the GLS sum 
rule; however, for other physical quantities measured at moderate energies the experimental 
errors may be much smaller, while the shift in central value would remain the same. 

We see therefore that within the modified perturbative approach it could be easier to 
accommodate in a consistent way some measurements giving relatively small value of the 
strong coupling parameter at moderate energies, as compared with the measurements of the 
strong coupling parameter at high energies. It would be interesting for example to analyze 
in detail within the modified approach the QCD corrections to the decay rates of the heavy 
quarkonia, which tend to give rather low values of the strong coupling constant, as has been 
recently discussed for example in |17C]. 



7. Summary and outlook 

In this article we considered the problem of reliability of perturbative QCD predictions 
at moderate values of (i.e. of the order of few GeV^). Using as an example the next- 
to-leading and next-to-next-to-leading expressions for the effective charge 6y, related to 
the static interquark potential, and the QCD correction <5gls to the Gross-Llewellyn- Smith 
sum rule in the deep inelastic neutrino-nucleon scattering, we have shown that QCD pre- 
dictions at moderate energies obtained with conventional renormalization group improved 
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perturbation expansion depend strongly on the choice of the renormahzation scheme (Sub- 
sections 2.3-2.5), even for a conservative choice of the renormahzation scheme parameters 
(condition ( ^.13D ). This casts doubt on the rehabihty of the commonly used perturbative 
expressions in this energy range. Taking closer look at the conventional expansion we ob- 
served that one of the possible sources of the strong renormahzation scheme dependence 
of perturbative predictions may be the strong scheme dependence of the conventionally 
used running QCD coupling parameter (the couplant) itself; an extreme manifestation of 
the scheme dependence of the conventional couplant is the presence of a scheme dependent 
singularity (Landau singularity) at nonzero value of (Subsection 3.1). 

As a step towards improving the reliability of perturbative predictions we proposed 
to replace the conventional couplant in the perturbative approximant of a given order by 
the modified couplant, obtained by integrating the renormahzation group equation with 
a modified, nonpolynomial ;9-function that in each order is a generalization of the poly- 
nomial expression for the /3-function used in the conventional expansion. We formulated 
some general criteria that the sequence of modified /3-functions has to satisfy in order to 
generate a useful modification of the conventional perturbation expansion (Subsection 3.2). 
We constructed a concrete model sequence of the modified /3-functions that satisfies our 
criteria, using the so called mapping method (Equation ( ^.13|) ). One of the properties of 
our model is that in each order of perturbation expansion the modified couplant is free 
from the Landau singularity, despite the fact that the /3-function does not contain any 
essentially nonperturbative terms. We then performed a detailed analysis of the modified 
NL and NNL order predictions for 6y and (5gls in various renormahzation schemes. In 
particular, we generalized to the case of the modified perturbation expansion the equations 
defining the renormahzation scheme distinguished by the Principle of Minimal Sensitivity 
(Subsection 4.2). The sequence of modified /9-functions in our approach contains two pa- 
rameters: parameter k, related to the asymptotic behaviour of the modified couplant in the 
limit 0, and the parameter r/, which characterizes the value of the modified couplant, 

for which the nonpolynomial character of the modified /3-function becomes important. The 
presence of these parameters gives us the opportunity to improve the accuracy of low order 
perturbative approximants by using use some information from outside of the perturbation 
theory. We proposed to fix the parameter k in such a way that the modified couplant in 
every order has a l/Q^ behaviour in the limit — > 0, and to adjust the parameter rj so 
that the modified NNL order PMS predictions for would match certain phenomenologi- 
cal expression for this effective charge (Subsection 4.3). We have then shown that for the 
chosen values of the these parameters the modified perturbative predictions for 5y and (^gls 
are much less sensitive to the choice of the renormahzation scheme than the predictions ob- 
tained in the conventional expansion (Subsection 4.4). The observed pattern of the scheme 
dependence of the modified expansion indicates in particular that the PMS predictions in 
this expansion stand on a better footing than in the conventional expansion. 

The general conditions on the sequence of the modified /3-functions in our approach 
leave some freedom in choosing the concrete form of these functions. In order to see, 
how this freedom might affect the modified predictions in low orders of the expansion we 
considered an alternative model sequence of /3-functions (Equation (|5.7|) ) and we repeated 
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the whole analysis of Sy and (Jqls for this new model (Subsection 5.3). We found that the 
results agree quite well with the predictions obtained in the previous model. We argued 
that in the PMS scheme this is likely to be true for any other /3-function satisfying our 
criteria. 

The modified predictions for 5^ and (5gls considered in this article have an interesting 
property that their evolution with is less rapid than that of the conventional predictions. 
In the case of the NL and NNL order PMS predictions we have given a general argument, 
that this would be true for all physical quantities in QCD of the form considered in this 
article, at least for relatively large values of (i.e. small values of the modified couplant). 
This has been achieved by expanding the PMS approximants in terms of certain regular 
couplant to one order beyond the formal order of the approximant (Subsection 6.1). To 
see, what might be the phenomenological significance of this effect we made a simple fit to 
the experimental data for 5gls (Subsection 6.2). After extrapolation to the energy scale 
of Mz we found that the obtained values of as{M'^) are shifted upwards (by about 5%) 
and are closer the world average (dominated by measurements at higher energies) than the 
values obtained with the conventional expansion. This suggests that within the modified 
expansion it might be easier to accommodate in a consistent way the relatively low values 
of the couplant obtained from some low-Q^ experiments and relatively large values of the 
couplant from some experiments at high Q^. 

In our discussion we have given explicit formulas describing the dependence of the 
modified couplant for various model /3-functions that are perturbatively consistent with 



the conventional expansion, but are free from Landau singularity (Equations ( 3.16 ), ( ^.91) , 
( |6.1| )). These very simple expression are interesting in their own right and might be useful 
in other theoretical considerations, unrelated to the concepts presented in this article. 

Our considerations in this article had exploratory character — they should be regarded 
more as a "feasibility study" of certain idea for the modification of the QCD perturba- 
tion expansion, than a very precise quantitative analysis of concrete predictions. It seems 
however that the results are sufficiently encouraging to justify a comprehensive study of 
various perturbative QCD predictions at moderate energies within the modified perturba- 
tion expansion approach. Such study could include the perturbative QCD effects in the 
physics of heavy quarkonia, in hadronic decays of the r leptons and in the deep inelastic 
lepton-nucleon scattering at the lower end of the energy scale. It would be then interesting 
to compare the results of such comprehensive analysis with the high energy measurements 
of the strong coupling parameter a^. As a first step towards such an analysis one should 
perform a more detailed study of possible methods to constrain the parameters rj and k in 
the proposed sequence of the modified /3-functions, generalizing the very simple approach 
adopted in this article. Another interesting line of research would be to study the resum- 
mation of the QCD perturbation expansion in the modified approach — the stability of 
the modified couplant with respect to change of the renormalization scheme should have 
favorable effect on the stability and the rate of convergence of the results obtained with 
various summation methods. In this context it would be interesting to investigate, whether 
the PMS prescription applied to the modified expansion does indeed act as an automatic 
resummation procedure for divergent perturbation series. 
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A. Parameterization of the renormalization scheme dependence of per- 
turbative approximants 



In this Appendix we give some more details on the renormalization scheme dependence of 
finite order perturbative QCD predictions. Let us begin with a non-RG-improved expression 
for the physical quantity 5 in the A^-th order of the expansion: 



+ f2{^i^Q^)a^N){^i^) + ...+ fN{^i^Q^)af,^){^^^) , (a.i) 



where // is the undetermined scale appearing in the process of renormalization and d(jv)(/"^) = 
is the couplant at the scale fi. The couplant satisfies satisfying the A^-th order 
renormalization group equation: 

= P^'^Ha.N)), (A.2) 

where /3^^^(d(^)) has the same form as in Equation (|2.2|). We introduced here the "hat" 
notation in order to distinguish the parameters fi{^'^, Q^) and d(/x^) from the corresponding 
parameters rj and a{Q'^) in the renormalization group improved expression for 5. 

The differential equation ([A.2|) is equivalent to the following transcendental equation 



^2 h \ 

-In-^ = Clin - + hcilno(m + F^^\a(]s,),C2, ...,cn), (A. 3) 

2 2 a(7v) 

where the function F^^\a,C2, ■■■,cisi) is given by ( |2.5|) . The parameter A appearing in 
Equation (|A.3| ) and the arbitrary integration constant have been chosen according to the 
conventions set by |^. 

The form of the coefficients fj is constrained by the fact that physical predictions of 
the theory should be (at least formally) independent of /i, i.e. 

M^^ = 0(a^^^). (A.4) 



Under our assumptions the coefficients fj may be written in the form 

n(/.^Q2) = r^°) + ^n^, (A.5) 



Q') = rf + (ci + 2rf^) ^ In ^ + Q In ^ V > etc. (A.6) 



where the parameters r^l^"^ are ordinary numbers. 

The couplant dx(Ai^) in some general renormalization scheme X is related to the MS 
couplant by a finite renormalization: 

aMs(^^) = «x(/i^) [1 + Aia^ii?) + A2al{i?) + ...] , (A.7) 
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where the constants Ai are related to the finite parts of the renormahzation constants and 
in principle may take arbitrary values. If the physical quantity 5 is expanded in terms of 
ax, the expansion coefficients take the form: 



-MS 



+ ^1, 



+ 2Aiff^^ + A2, etc. 



(A.8) 



In general in the new scheme also the coefficients in the /3-function would be different; for 
example in the NNL order we have 



C2 



cf^ + A,c, + Aj 



The change of the renormahzation scheme affects also the parameter A ||5|]: 

Ax = Aj;jgexp 



(A.9) 



(A.IO) 



This relation is exact to all orders of the perturbation expansion. 

In the A^-th order of the perturbation expansion we have A^ parameters (Ai) character- 
izing the freedom of choice of the renormahzation scheme and 2N — 1 expansion coefficients 
(rj and c^), which depend on these parameters. It is clear therefore that the expansion 
coefficients are not totally independent — for N >2 there must be A^ — 1 combinations of 
these coefficients, which are independent of the renormahzation scheme |jll|, ^, |l3|, [ij, [T5| . 
For N = 2 the relevant scheme invariant combination may be written in the form: 



P2 = C2 + f2 - Cifi - fi . 



(A.ll) 



It should be noted that the invariant shown above corresponds to the definition adopted 
in 1^, |l^, U, |l5|], which is slightly different from the definition introduced in |11|. It is not 
surprising that we have some freedom in defining the scheme invariant combinations of the 
expansion parameters, because any combination of the scheme invariants is of course itself 
a scheme invariant. It seems, however, that the invariants defined by P, 13, 14, ll^ give a 



more natural measure of the magnitude of the higher order radiative corrections. This may 
be seen in the following way: Let us take for example the NNL order approximant S^'^\ 
calculate the derivative Q"^ d5^'^\Q'^) / dQ'^ and expand the result in terms of 5^^-* itself. We 
obtain 



k=3 



(A.12) 



This looks very much like the renormahzation group equation (^), except that instead of 
the coefficient C2 we have the scheme invariant combination p2 (exactly in the form proposed 
in 1^, 13, 14, |l^), while the higher order expansion coefficients pk in this expression are — 
of course — scheme dependent. If we would do the same calculation with 5^^\ we would 
obtain a similar equation, where the coefficients up to and including the order N = 3 would 
be the scheme invariants pi defined by |||, |l3|, U, |l5|, while the coefficients of higher order 
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would be scheme dependent. This shows that the invariants proposed in H, |T^, U, ^ have 
indeed some universal meaning. 

In phenomenological applications it is usually assumed that renormalization scale fi is 
proportional to the characteristic energy scale of the process: fjp = A^Q^, where A is some 
constant. In this way we obtain the renormalization group improved expression ( p] ) for 
the physical quantity 6, where the coefficients rj are now independent of and in the 
arbitrary scheme X they take the form 

n = rf^^ + blnX + Ai, (A.13) 
r, = ^ + (ci + 2ri°)^) 61n A + (61n A)^ + 

+ 2Ai{rf'^^^ + blnX)+A2, etc. (A. 14) 



The whole Q^-dependence of 6 comes then from the Q^-dependence of the couplant a(^j^-^{Q'^ 



a/]^\{}?Q'^), which is determined by the implicit equation 



5 q2 ^ 
- In — = -Ai - 6 In A + ci In - + 

2 2 



+ + ci lna(7v) + F^^^(a(7v),C2, ...,C7v), (A. 15) 



where we used the exact relation ( A. 10 ) to introduce A^^ as a reference phenomenological 



parameter. The use of a{Q'^) instead of a{iJ?') in ( |2.1D is equivalent to resummation of some 
of the (In^)*^ terms appearing in rfc(/i^,Q^) to all orders. Of course a(jv)(Q^) satisfies 
Equation ([2^). 

To obtain a convenient parameterization of the available degrees of freedom in choosing 
the perturbative approximants we first note, that the next-to-leading order approximants 
contain in principle two arbitrary parameters, Ai and A, which however appear in the 
expression for 5^^^ and in the Equation ( ^.41 ) in the combination Ai + 6 In A. This means 
that in the NL order the freedom of choice of the approximants may be characterized by 
only one parameter; we found it convenient to choose as such a parameter the coefficient 
ri in the renormalization group improved expression for 5. Thus in the NL order we have: 

5^^\Q\n) = a{Q\n) [l + na{Q\ri)\ , (A.16) 

with a(Q^,ri) determined by the equation (^) with 

F(^)(o) = -ciln(l + cia). (A.17) 

It should be emphasized, however, that although under our assumptions the parameters 
Ai and A have the same effect on 5^^^ , they have in fact a completely different character. For 
example, in a more general class of renormalization schemes the coefficient Ai may depend 
on quark masses and the gauge parameter. The procedure of fixing the choice of the 
renormalization scheme (i.e. Ai) differs from the procedure of fixing the renormalization 
scale (A), even in the NL order. This observation has important consequences for the 
discussion of the "reasonable" schemes for NL order approximants. For example, it shows 
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that one should be very careful discussing "reasonable" value of the scale parameter; as 
has been emphasized for example in |17|, a given value of ri — which provides a unique 
specification of the NL order approximant — may correspond to a "reasonable" value of A 
with one subtraction procedure (say, MS), and a completely "unreasonable" value of A in a 
scheme defined by some other subtraction procedure (for example momentum subtraction). 

In the next-to-next-to-leading (NNL) order there appears an additional degree of free- 
dom in choosing the renormalization scheme, corresponding the freedom of choice of the 



parameter ^2- Following ||l]J we shall parameterize this degree of freedom by the coeffi- 
cient C2 in the NNL order /3-function. If the parameters ri and C2 are fixed, the value of 
the coefficient r2 may be determined from the scheme invariant combination p2: 



r2(r-i, C2) = P2-C2+ cm + rl. 



(A.18) 



Thus in the NNL order we have: 

<5(2)(Q2^ri,C2) =a(Q2,ri,C2) [l + ria(Q^ n, C2) + r2{ri,C2) a^{Q\n,C2)\ , (A.19) 



where a{Q'^ ,ri,C2) is determined by Equation (2^) with F^'^\a,C2) of the form (for C2 > 
cf/4) 



F(2)(a,c2) = -^ln(l + cia+" 



C2a ) + 



2C2 



arctan 



2 + cia 



(A.20) 



The expression for F^'^\a, C2) for other values of C2 may obtained via analytic continuation 
in C2. The coefficients of the /^-function have the following values: b = (33 — 2nj)/6 |[l|, ^, 
Cl = (153 - 19n/)/(66 - 4n/) |7l|, |172|, [17| and [[7|, [17| 



MS 
C2 



77139 - 15099 n/ + 325 n2 



288(33 - 2nf) 



B. Principle of Minimal Sensitivity scheme for conventional expansion 

In this Appendix we discuss explicit formulas for the PMS parameters in the conventional 
expansion. It should be noted, that the parameterization of the RS dependence and the 
form of the NNL order scheme invariant p2 adopted in this article are different from those 



assumed in the original paper |11|. 

In the NL order the scheme parameter n corresponding to the PMS scheme is a solution 



of the Equation (2^. Performing the differentiation and taking into account that in the NL 
order we have 



(B.l) 



we obtain 



a2 + (l + 2fia) -/?«(( 



Solving this equation for ri we find 



Cl 



2(1 +cia)' 



(B.2) 



(B.3) 
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Inserting this expression into the Equation (|2.4| ) for the NL order couplant we obtain a 
transcendental equation for a. Solving this equation we obtain the numerical value of a for 
the chosen value of , which then allows us to determine the numerical value of fi , and 
hence the numerical value of 6^^^ in the PMS scheme for this Q^. For small a (i.e. for large 
(5^) we may use an approximate expression for fi: 

ri = -y + 0(a). (B.4) 

In the NNL order the parameters fi and C2 corresponding to the PMS scheme are 
solutions of the system of two equations (|2.7|) and (|2.8|). In the NNL order we have 



so the partial derivative over ri has the form: 

AjCs) = ^2 ^ ^ 2n)a^ + (1 + 2na + 3r2a2)|/3(2) („) (b.6) 

= -(2ciri + C2 + 3r2)a^ - (2riC2 + 3cir2)a^ - 3r2C2a^ (B.7) 

The equation ( |2.7D is therefore equivalent to: 

3p2 - 2c2 + 5cifi + 3rl + (2fiC2 + 3f2Ci)a + 3f2C2a^ = 0, (B.8) 

where in the 0{aP) term we expressed f2 in terms of n and C2. The partial derivative of 
over C2 has the form: 

^^(2) = _„3 + (1 + 2ria + 3r2a2)|^ , (B.9) 

OC2 OC2 



where 



More explicitly 



dc2 b'^ ' dc2 Jo (/3(2))2 dc2 



dc2 

where for C2 > ^ we have 

4C2 



— = a^ + H{a,C2), (B.ll) 



H{a,C2) 



(4C2 - c2)3/2 



2/1 2n a a/ 4c2 

a (1 + cio + C2a j arctan 



2 + cia 



2 ^ 2 ^ 



(B.12) 



For other values of C2 the relevant expression is obtained via analytic continuation in C2. It 
is easy to verify that H{a,C2) = 0{a^). The equation (|2.8|) for the PMS parameters may 
be therefore written in the form 

2fi + 3f2a + (1 + 2fia + 3f2a2) = 0. (B.13) 
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Solving the PMS equations ([B.8|) and ( |B.13|) we obtain the parameters fi and C2 singled 
out by the PMS method, expressed in terms of a. Inserting this solution into the implicit 
equation for the NNL order couplant ( |2.4|) we obtain the parameter a for the chosen value 
of . Inserting these values into the expression ( |A.19 ) we then obtain the PMS prediction 
for (5(2) at this Q^. 

The equations ( |B.8| ) and ( [B.13| ) are quite complicated, but for small o (i.e. for large 
Q2) it is easy to solve them in an approximate way ||176| , |177| , |178| , ITS]. If we look for fi 
and C2 in the form of a series expansion in o, then from ([B.13| ) we see, that 

fi = 0{a). (B.14) 

From ([B.8| ) we then immediately find that 

C2 = \p2 + 0{a). (B.15) 
C. Solution of the inequality cr2(ri,C2) < l\p2\ 

In this Appendix we give a detailed description of the sets of the NNL order renormalization 
scheme parameters ri and C2) that satisfy the condition ( 2.13| ) 

a2{ri,C2)<l\p2l (C.l) 

where / > 1. Let us recall, that p2 is given by the Equation ( |2.3D 

P2=C2+r2- cm - rf, 

and a2{ri,C2) is given by the Equation ( p. 12 ) 

o-2(n,C2) = |c2| + |r2(ri,C2)| + ci|ri| +rl. 

We assume of course p2 / 0. We have to consider several cases, depending on the value of 
P2 and I. 

1. For p2 > cf/4 we first define the following parameters: 



rr = -Vp2(/-1)/2, (C.2) 

= I {-ci + ^cf + 2(/-l)p2) , (C.3) 

cf- = -p,{l - l)/2, (C.4) 

c-- = p2(/ + l)/2, (C.5) 

cf = cirf" + cr". (C.6) 

For C2 > the set of scheme parameters in the (ri, C2) plane satisfying the condition ( \C.l\ ) 

is bounded by the segments joining points (r™™,0), (r™™,4°*), (0,0^"^), {r"^"^^ , df^^) , 
(r™^^,0). For C2 < this set is bounded by the curves 

C2(ri) = rl + cf" for rf" < n < 0, (C.7) 
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c^^n) = rl + cm + cf" for < n < rP"- (C.8) 

2. In the case < < cf/4 and I ^< 4-rf-4v^ ^ 4-^+4^1=3 ^j^^^.^ ^ ^ 4p2/cf , the 
description of the relevant set of parameters is the same as in the case p2 > cj/A. For 
^ i-d-iVT^ ^ > have to redefine the parameters rf"" and 4^^*: 

rf- = -p,{l + l)/2ci (C.9) 
c^' = {rf'f + 0^". (C.IO) 



For C2 > the set of parameters satisfying the condition (|C.lj ) is bounded by the segments 
joining points (rf™,0), (OjCf"^), (rf c™'^''), For C2 < this set is bounded by 

the curves 

C2(ri) = rj + cf" for rj""^ < n < 0, (C.ll) 
C2{ri) = rl + cm + cf"" for < n < rr^", (C.12) 
supplemented by a segment joining points (r™™,0) and (r™™,c™*). 

3. For P2 < and I/O2I > 2cf(/ + !)/(/ — 1)^ we define the parameters 



^1 



-V|p2|(/ + l)/2, (C.13) 
1 



2 1 -ci + Vcf + 2a + l)|p2|), (C.14) 

cr = -|P2|(/ + l)/2, (C.15) 

= |P2|(/-l)/2, (C.16) 

= cirj"'"^ + cr"- (C.17) 



For C2 > the set of scheme parameters in the (ri, C2) plane satisfying the condition ( C.ip 
is bounded by the segments joining points (r™™,0), (r™™,C2°*), (0, c™'^''), (r™^"", c™^''), 
^^max^g-j p^j, C2 < this Set is bounded by the curves 

C2(ri) = rl + cf^'"^ for r™'^ < ri < 0, (CIS) 

C2{ri) = rl + cm + cf'' for < n < rr^"- (C.19) 

4. For p2 < and \p2\ < '2.c\{l + !)/(/ — 1)^ we have to redefine again some parameters: 

r-- = -|p2|(/-l)/2ci, (C.20) 

^int ^ (^min)2 ^ ^min^ 21) 



For C2 > the set of scheme parameters in the (ri,C2) plane satisfying the condition (|C.1| ) 
is bounded by the segments joining points (r™™,0), (0, c™'^''), {r^^^.c^^^), (r^'^'^jO). For 
C2 < this set is bounded by the segment joining points (r™°,0) and (r™™,C2°*) and the 
curves 

C2(ri) = rl + cf'"^ for rf'' < ri < 0, (C.22) 
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c^^n) = rf + cm + cf" for < n < rP". (C.23) 

One way to obtain these solutions is to consider a three dimensional space (ri,C2,r2) 
and carefully analyze all possible intersections of a (closed) surface defined by the equation 

|C2| + |r2| +ci|ri| +rj = 2|/92| (C.24) 

with the surface 

1^2 = P2 - C2 + ciri + rf. (C.25) 
Calculations are straightforward, but tedious. 
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